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ABSTRACT 

A  pneumatic,  pure  fluid  oontrol  system,  capable  of  main- 
taining the  angular  position  of  an  inertia  load  in  a  gravita- 
tional space,  using  only  one  moving  part,  was  designed  and 
built.  The  control  system  consisted  of  a  vortex-sink  angular 
rate  sensor,  a  pendulum  controlled  flapper-nozzle  angular  posi- 
tion sensor,  a  summing  amplifier  and  bistable  reaction  jets. 
The  system  was  designed  to  operate  in  a  bang- bang  manner  in 
order  to  reduce  errors  in  a  minimum  amount  of  time. 

A  previous  author's  work  on  the  vortex-sink  angular  rate 
sensor  is  extended  to  include  information  and  data  on  propor- 
tional amplification  of  the  output  signal  and  the  frequency 
response  of  the  sensor.  Methods  for  synthesizing  the  components 
into  a  working  system  are  presented  along  with  a  theoretical 
analysis  of  the  system's  performance. 

Results  of  experiments  performed  on  the  completed  system 
showed  that  the  reaotion  jets  could  be  switched  in  a  nearly  op- 
timal manner  when  responding  to  errors  introduced  on  one  side 
of  the  system,  but,  due  to  the  unsymmetrical  gain  of  the  rate 
sensor,  the  errors  introduced  on  the  opposite  side  resulted  in 
a  certain  amount  of  after  end-point  chattering. 

Good  correlation  between  theory  and  experiment  led  to  the 
conclusion  that  the  application  of  standard  control  system  design 
techniques  oan  be  readily  applied  to  the  design  of  pure  fluid  oon- 
trol systems. 
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CHAPTER  I 
INTRODUCTION 

1.1  General  Introduction  to  Pure  Fluid  State  Technology 

The  basic  concept  of  pure  fluid  amplification  is  that 
high  energy  fluid  streams  can  be  controlled  by  low  energy 
streams.  The  concept  was  evolved  at  the  Harry  Diamond  Lab- 
oratories of  the  U.  3.  Army  Material  Command,  "Washington, 
D.  C.  in  1959.      ' 

At  the  First  Fluid  Amplification  Sympqsium  held  at 
the  Laboratories  in  October  of  1962,  Brigadier  General 
John  G.  Zierdt  described  the  possible  future  of  these  compon- 
ents : 

This  advance  has  removed  many  of 
the  limitations  on  conventional 
military  end  commercial  pneumatic 
systems  and  has  enabled  an  extreme- 
ly rapid  extension  of  their  capabil- 
ities.  The  use  of  low  energy  flow 
to  direct  and  control  high  energy 
streams  permits  the  design  of  pure 
pneumatic  components,  elements,  and 
systems,  capable  of  amplification, 
feedback,  memory,  logic  functions, 
analogue  computation,  and  digital 
computation.  This,  in  turn,  makes 
it  possible  to  use,,  pneumatic  control 
systems,  amplifiers,  and  computers 
where  they  could  not  be  used  before. 

1.2  Purpose  and  Scope  of  Thesis  , 

Sinoe  the  inception  of  pure  fluid  amplification,  much 


work  has  been  done  in  analyzing  individual  components,   re- 
cently the  emphasis  in  the  field  has  been  turned  to  the  de- 
velopment of  complete  systems.   Although  much  work  has  been 
done  along  these  lines,  the  quantitative  data  associated  with 
such  systems  is  scanty  and  inaccessible. 

One  objective  of  this  thesis  is  to  develop  techniques 
that  may  b«.  used  in  designing,  analyzing  and  synthesizing  these 
components,  into  a  working  control  system. 

The  specific  objective  is  to  develop  a  pneumatic  con- 
trol system,  using  as  few  moving  parts  as  possible,  capable 
of  maintaining  a  commanded  angular  position  of  a  pure  inertia 
load.  Should  errors  be  introduced,  the  control  system  is  to 
reduce  these  errors  to  near  zero  in  the  shortest  possible 
time  within  the  limitations  of  the  prime  move  used. 

Another  design  objective  is  to  purposely  make  the 
system  respond  slowly  so  that  the  device  may  be  used  as  an 
educational  model  demonstrating  the  working  of  individual 
components  and  the  system  as  a  whole. 

1.3  Discussion  of  the  Proposed  System  to  Attain  the  Design 
Objectives 

Minimum  time  control  of  a  vehicle  generally  involves 
the  use  of  th£  full  force  available  from  the  prime  mover  in 
a  bang-bang  manner,  jfor  example,  it  is  obvious  that  to  move 
an  automobile  from  a  point,  A,  to  another  point,  B,  in  a  min- 
imum amount  of  time,  the  operator  should  use  maximum  acceler- 
ation, at  the  right  instant  lock  his  brakes  and  slide  into  the 


desired  end  point.  The  operator  in  this  example  would  be 
using  full  driving  force  in  one  direction  and  full  braking 
force  in  the  other  to  achieve  the  minimum  time  objective. 
In  controlling  an  angular  inertia  load,  it  con  be 

shown t that  minimum  time  movements  are  obtained  by  using  full 

1  * 
torque  in  one  direction  or  the  other  .  l^gure  1.1  shows  a 

phase  plane  portrait  of  an  inertia  load  being  driven  by  a 

constant  torque,  whose  sense  is  either  plus  or  minus,  for 

various  initial  conditions. 

It  can  be  seen  from  Figure  1.1  that  one  curve  exists, 
which,  if  used  as  a  boundary  to  switch  the  sense  of  the  torque^ 
will  allow  the  system  to  reduce  any  errors  to  zero  in  mini- 
mum time.  In  order  to  switch  the  sense  of  the  applied*  torque 
at  the  proper  instant,  it  is  obvifcus  that  the  control  system 
must  be  able  to  l)  measure  angular  positions,  2 J  measure 
angular  rate,  3)  compute  the  coordinates  of  the  switching 
curve,  and  4 J  compare  the  values  of  angular  rate  and  angular 
position  with  the  values  of  the  switching  curve  and  switch 
the  sense  of  the  applied  torque  at  the  proper  instant. 

The  remaining  chapters  of  this  thesis  deal  with  the 
design  of  the  components  needed  to  accomplish  the  above  men- 
tioned functions  and  the  synthesis  of  these  components  into 
a  working  system.  , 
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Figure  1.1 

Phase  Plane  Representation  of  an  Inertia  Load,  being 

Driven  by  a  Constant » Torque 


CHAPTER  II 
VORTEX  ^NGUL^.R  RATE  SENSOR  ANALYSIS  AND  DESIGN 

2.1  General  Theory  of  Operation 

The  principles  on  which  a  vortex  angular  rate  sensor 
operate,  are  quite  simple.   Real  difficulties  quickly  arise, 
however,  when  attempts  are  made  at  mathematical  analysis  of 
the  complex  flow  which  takes  place  in  these  devices  or  when 
attempts  are  made  to  obtain  pressure  or  flow  signals  which  are 
some  function  of  the  input  angular  rotation.   It  is  primarily 
because  of  the  latter  difficulty  that  the  methods  of  obtain- 
ing output  signals  have  been,  and  still  are,  held  confidential, 

The  basic  operation  of  a  rate  sensor  is  depicted  in 
figures  2.]aend  2.1b.   Figure  2.1a  shows  the  flow  state  in  the 
absence  of  angular  rotation.  In  this  case  a  particle  of 
fluid  passes  from  the  supply  plenum  through  the  porous  ring 
into  the  vortex  chamber  and  out  the  exit  hole.  From  mass 
continuity  the  radial  flow  velocity  is  described  by 

Ur   -    — 9. 12.1) 

2-rrr  h 

where 

^  ■  flow  rate  supplied  to  the  sensor 

r  =  any  arbitrary  radius 

h  =  thickness  of  the  vortex  chamber 
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;ure  2.1 
Vortex  liar    la  to  3  en.  s  or 


Thus  just  inside  the  porous  ring,  the  radial  velocity 
is  given  by 


and  at  the  exit 


Ue  =  cl (2.3) 

2-«Tiieh 


where  Q,  and  h  are  as  given  above 
and  HQ  =  radius  of  porous  ring 
H0  =  exit  hole  radius 
Figure  2.1b  shows  the  flow  condition  when  the  sensor 
is  rotated  at  a  constant  angular  rate.  In  this  case  a  par- 
ticle of  fluid  entering  the  vortex  chamber  has  the  radial 
velocity  given  by  equation  (2.1)  and  a  tangential  component 
of  velocity  imparted  to  it  while  passing  through  the  porous 
ring.  In  visualizing  this  tangential  component,  one  must 
use  a  frame  of  reference  other  than  the  vortex  chamber.  Thus 
with  respect  to  laboratory  coordinates,  the  magnitude  of  the 
tangential  component  of  velocity  is  given  by 

Ut0  =  a0GJ  (2.4) 

where  U-to  =  magnitude  of  the  tangential  velocity 
CJ  ■  magnitude  of  the  input  rotation 
Hence,  a  particle  of  fluid  of  mass  dm  will  have  angular  mo- 
mentum given  by  3q  ^to  dm*   Conservation  of  Angular  Momentum 
produces  an  expression  for  the  tangential  velocity  (Ut)  at 
an-  arbitrary  radius  (r) : 

r  Ut  dm  ■  R0  Uto  &a.  (2.5a) 

Ut  »  Ro2a)  (2.5b) 


For  an  ideal  flow  the  familiar  result  of  Ut  approaching  infin- 
ity is  seen  as  r  approaches  zero. 

Thus  it  can  be  seen  that  if  a  means  can  be  devised  to  meas- 
ure the  amount  of  tangential  velocity  directly  or  indirectly  at 
a  particular  radius,  this  velocity  will  change  in  proportion  to 
the  magnitude  of  U)  as  to  the  direction  of  rotation. 

The  method  used  to  measure^  the  amount  of  tangential  veloc- 
ity present  in  the  vortex  chamber  of  the  rate  sensor  used  in  this 
control  system  was  developed  by  Sarpkaya4.   The  analytical  analy- 
sis about  to  be  discussed  is  also  given  by  Sarpkaya  but  is  included 
here  to  maintain  continuity  of  this  chapter  and  to  set  the  back- 
ground for  the  extension  of  Sarpkaya *s  work  which  follows  in  the 
remaining  articles*.  , 

A  schematic  of  the  sensor  developed  by  Sarpkaya  is  shown 
in  Figure  2.2.   This  sensor  differs  from  the  one  shown  in  Figure 
2.1  in  that  two  sink-tubes  are  attached  to  the  exit  holes  of  the 
vortex  chamber.   These  sink-tubes  serve  the  following  purposes: 

1)  Tangential  flow  amplification  -  each  tube  can  now  be 
considered  as  a  miniature  elongated  vortex  chamber  being  fed  from 
the  top  by  a  uniform  radial  flow  on  which  is  superimposed  the 
tangential  swird  flow  from  the  main  vortex  chamber, 

2)  Directing  radial  flow  -  the. tubes  receive  the  flow 
from  the  main  vortex  chamber  and  direct  it  toward  the  probe, 

3)  The  sink-tubes  provide  a  housing  for  the  probe. 

The  probe,  shown  in  Figure  2.3,  consists  of  a  small  cir- 
cular cylinder  divided  into  two  chambers.  Two  small  holes  are 
drilled  perpendicular  to  the  center  line  of  the  probe  and 


*  To  further  credit  the  work  of  Sarpkaya,  most  of  his 
symbols  are  maintained  throughout  this  chapter. 
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Figure  2.3 
Cross  Section  of  Probe 
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aligned  along  a  straight  line  drawn  on  the  surface  of  the 
cylinder  between  its  ends.   The  probe  is  inserted  into  the 
sink-tube  perpendicular  to  the  axis  of  flow,  as  shown  in 
Figure  2.2.   The  holes  are  set  at  an  angle  to  the  oncoming 
flow  of  fluid  in  the  sink-tube. 

The  operation  of  this  probe  in  sensing  angular  rota- 
tion is  as  follows:   Jfrom  potential  flow  theory  the  pressure 
distribution  around  a' circular  cylinder  placed  in  a  uniform 
flow  field  is  given  by 

?b  '     Po  +  PU5     (1  "  4sin2$     (2.6) 
a 

where  Pj>  s     pressure  in  lbf/in2 

Po  =  static  pressure  in  lbf/in2 

O  -  density  in  lbf-sec2 

in4 

Us  ■  axial  sink  velocity  in  in/sec 

(/    =     angular  position  measured  eround  the 
cylinder  from  the  stagnation  point . 
A  sketch  of  equation  (2.6)  is  shown  in  figure  2.4.   in  the 
absence  of  angular*  rotation  of  the  device,  the  uniform  radial 
velocity  in  the  vortex  chamber  (Ur)  turns  into  the  sink-tube 
(now  redesignated  Us) ,  strikes  the  probes  and  then  passes  out 
of  the  sensor.   Because  both  holes  are  aligned  at  the  same 
angular  position  along  the  probe,  equal  pressures  exist  in 
each  chamber.  The  output  signal  of  the  probe,  however,  is 
the  difference  between  these  two  pressures.   Sinoe  both 
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Pressure  in  probe 
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ideal  angle  of  incidence 
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Figrure  2.lf 


Ideal   Pressure  Distribution   Ground  a  Circular  Cylinder 
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chambers  are  at  the  same  pressure,  no  output  signal  (Ai") 
exists. 

Nov;  let  there  exist  angular  rotation  of  the  rate  sensor. 
Figure  2.5  shows  the  velocity  vectors  in  way  of  the  probe. 
The  tangential  velocity  (Vp0)  now  present  in  the  sink-tube 
has  the  feffeot  at  diverting  the  main  axial  sink  flow.  The 
top  hole  is  thus  placed  in  a  uniform  flow  coming  from  an 
angle  which  is  greater  than  that  which  existed  without  angu- 
lar rotation.  The  bottom  hole  is  placed  in  a  flow  coming  from 
an  angle  which  is  smaller.  Thus,  the  pressure  in  the  chamber, 
to  which  the  top  hole  is  connected,  decreases, (point  A  in 
Figure  2.4),  while  the  pressure  in  the  bottom  chamber  increases, 
(point  B  in  Figure  2.4),  A  pressure  differential  now  exists 
across  the  probe  which  is  a  function  of  the  magnitude  of  the 
angular  rotation.  It  is  readily  visualized  that  had  the  angu- 
lar rotation  been  in  the  opposite  direction  that  the  reverse 
procedure  would  have  taken  place.  Hence,  the  signal  also 
senses  direction  of  angular  rotation  as  well  as  magnitude. 
It  is  assumed  in  this  explanation  that  the  relative  magnitude 
of  the  sink-tube  axial  and  tangential  velocities  are  such 
that  a  long  spiraling  flow  exists  at  the  radius  rpo. 

To  visualize  the  flow  quantitatively  and  to  obtain  the 
function  relating  Ap  an^  COrefer  again  to  Figure  2.5.   By 
.inspection  it  can  be  seen  that  A$,  Us,  and  vpo  are  related 

A#=  tan"1   VPQ  (2.7a) 

Us 
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sink- tube 


probe 


Figure  2.5 
Sketch  of  Probe  in  Operation 
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If  A$  is   a   small  angle   then: 


A#      =  ^  (2.7b) 

,    To  obtain  the  cazimum  signal  (A.P)  possible  for  a 
given  Cc),it  is  necessary  to  determine  the  maximum  slope  of 
the  curve  shown  in  Figure  2.4.  Hence,  cif ferentiating 
equation  (2.6)  to  obtain  the  slope  ^ives: 

^Pb         Pus2   SSin^Cos^  (2.8) 

30    "   "    2 
Now  the  derivative  of  equation  (2.8)  set  to  zero  will  give 

the  0  for  maximum  slope. 


2 

e 

30  *~  2~ 


2  A  T 

3  ft   =    -  8   PUs    cos  2$,^tr/4 

1   *c 


(2.9) 

Therefore,  the  probe  holes  should  be  located  at  a 
45  degree  angle  to  the  oncoming  flow.   This  value  substituted 
into  equation  (2.8)  gives 

a  Pb     =    -  4  /pu32   A  (2.10) 

and   in  differential  form  for  one  hole: 

API        =  4  f  PUs        \    A£  (2.11) 


Equation   (2.5b)    gives  the  relation  between  the   tangential 
velocity  existing  in  the  sink-tube   and  the   input   angular 
rotation.      Y/ith  new  symbols,    equation    (2.5b)    becomes: 

V     =     *oaa>  (2-12> 

rpo 
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Combining  equations  (2.7b),  and  (2.11),  and  (2.12)  gives: 


A  p   1 

-"     4     Rn2W 

pUs2 

0 

Usrpo 

»  for  one  hole  (2.13a) 


Sinoe  two  holes  are   being  used  equation   (2.1X3a)    beoomes: 


Ap 


Pu 


=     8     BQ~  CO 


(2.13b) 


Usr  „ 
s    po 


This  is  the  desired  relationship  between  Ap  and  CO  in  non- 
dimensional  form.  Note  that  for  a  given  sensor  and  probe  con- 
figuration, that  R0  and  r   are  fixed.  Further,  if  the  flow 
through  the  sensor  is  fixed,  then  Us  is  constant.   Thus  it  is 
seen  that  equation  (2.13b)  is  of  the  form: 


Ap 


KOi 


($.14) 


where       X       =         8       pu  ; 


HQ6  CO       ,   a   constant 
—     depending  on 


2  ^srpo  geometry  and 

flow. 

Data  given  by  Sarplcaya  bear  out  the  linearity  of  this  equation 

over  a  substantial  range  of  CO   .  As  might  well  be  expected, 

from  inspection  of  Figure  2.4,  linearity  is  preserved  provided 

A  (7  is  small. 

Thus  far  the  equations  derived  to  describe  the  opera- 
tion of  this  rate  sensor,  are  based  on  an  idealized  potential 
flow.  It  is  now  appropriate  to  consider  the  effects  of  vis- 
cosity  on  the  operation  of  the  sensor  and  the  describing  equa- 
tions. Two  major  effects  of  viscosity  will  be  noticed: 

1)   Viscosity  will  cause  the  vortex  motion  in  the 
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vortex  chamber  and  sink -tube  to  deviate  from  the  ideal  vortex 
give:,  jy  equation  (2.5b). 

2)   Viscosity  will  cause  the  flow  around  the  probe  to 
separate  resulting  in  an  alteration  of  the  pressure  distri- 
but  ion  around  the  pro'oe  given  by  equation  (2.6). 
Effect  1)  is  studied  in  great  detail  by  Sarpkaya,  by  analyzing 
the  boundary  layer  equations  in  the  chamber  for  the  radial  flow 
with  small  tangential  flows  superimposed  on  the  main  flow. 
The  results  show  that  viscosity  reduces  the  strength  of  circula- 
tion as  the  flow  travels  from  R0  to  rp0.  Sarpkaya  introduces 
a  factor  he  calls  the  viscous  efficiency  of  the  sensor  defined 
as: 

E   =  -JiJE (2.15) 

To 

where  fp  =  strength  of  vortex  circulation  at 
radius  r 

P0  =  strength  of  vortex  circulation  at 
radius  R0 

The  viscous  efficiency  between  the  porous  ring  and  the  probe 

is: 

p 
E     =  P°  =     vdo   2^rPo        =     vpo   rPQ 

^°  ^o^2tTH0  *o  & 

(2.16a) 

from  which 

v        =  Ro   ^      E  (2.16b) 

rpo 
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Wj.th  equation  (2.16b)  used  in  lieu  of  equation  (2.12)  and  then 
combined  with  eou-ations  (2.7b)  and  (2.11),  the  new  relation- 
ship between  AP  and  CO  becomes 

IA4*    -  8  Hn2^E   for  two  holes.   (2.17) 
2 

# 

Intuitively  one  might  expect  thi.s  efficiency  to  be  relatively 
high  provided  the  circulation  is  quite  small.  This  is  shown  to 
be  true  by  Sarpkaya  in  a  plot  of  E  versus  Reynolds  Number. 

Effect  2)  is  an  extremely  important  consideration  which 
seems  to  have  been  neglected  by  Sarpkaya  in  his  analysis  of 
the  sensor.  Non  ideal  flow  around  the  probe  radically  re- 
duces the  slope  of  A  P  versus  Q  in  the  vicinity  of  $-  45°. 

For  the  case  of  ideal  flow  around  the  probe,  the  nondi- 
mensional  slope  is  given  by  equation  (2.11). 

=  ^  (2.11) 


a  /n  repeated 

From  actual  measurements  taken  on  the  probe  used  in 
the  sensor  constructed  by  the  author,  this  slope  was  found  to 
vary  between  2.58  and  3.17  for  the  flows  considered.   The  ef- 
fect of  viscosity  is  more  clearly  shown  in  Figure  2.6  where 
nondimensional  pressure1  distribution  for  two  actual  flows 
versus  (j     is  shown  along  with  that  predicted  by  potential 
flow  theory.  Thus  the  relationship  between  AP  and  CO  with 
the  introduction  of  the  new  slope  of  Ap  vs  A  (?   becomes: 
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—  TT0  =   171  .9   ft/sec 


U     =   199.5  ft/sec 
s 


&  (degrees) 


Figure  2.6 

Effect  of  Viscosity  on  the  Pressure  Distribution 

Around   a   Circular  Cylinder 
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AP  j    =   5.16  Ro2c°Z  for   tv;Q  h0le3# 


jOUs2/2         Us  rp0  (2.18) 

2.2  Rate  Sensor  Design  and  Design  Procedure* 

The  author  did  not  consider  the  effect  of  slope  reduc- 
tion due  to  viscosity  of  the  AP  versus  Q  curve,  discussed 
at  the  end  of  article  2.1,  until  after  the  rate  sensor  had 
been  designed.  Kence,  the  equations  used  for  the  design  were 

AP  =   elP^s2      j  &Qt    two  holes   (2.19) 

A$  =     V  CQE  (2.20) 

rpo  Us 

i 

From  the  design  specifications  discussed  in  the  intro- 
ductory chapter  to  this  thesis,   ^^iax  was  selected  to  be 
70  degrees  per  second.  In  order  that  the  output  of  the  sen- 
sor be  linear  over  the  range  -W^j-W  £ CO^qz*    tlie  vaIue  °^ 
A$  must  be  constrained  to  stay  within  the  linear  portion 
of  the  AP  versus  Q   curve.  Analysis  of  the  data  given  by 
Sarpkaya  has  shown  that  A  Q   should  not  exceed  roughly ±3° 
for  linearity. 

Inspection  of  equation  (2.20)  shows  that  RQ,  rp0,  E, 
and  Us  are  still  to  be  determined.  RQ  was  chosen  on  size 
considerations  of  the  sensor  and  was  selected  to  be  three 
inches.  The  value  of  Us  will  depend  on  the  magnitude  of  the 


*  See  Appendix  A.l  for  detailed  numerical  calculations 
in  this  section. 
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desired  output   signal  A?  and   can  be   determined  from     equa- 
tion    (2.19).     Once  Us   is   determined,   **  can  be   calculated  from 
the  relationship 


-[ 


ii/2     =       oTrg2  -  2c(2rs) 


Uc 


(2.21) 


where  r3  -    radius  of  sink 

c   =    radius  of  probe 

In  order  to  limit  Q,  to  a  reasonable  value,  it  is  obvious  that 
the  magnitude  of  A?  must  be  kept  relatively  small.   For  this 
design  a  value  of  Ap  »  0.1  psi  to  occur  at  CJS  ^max  was 
selected. 

Having  determined  ^,  the  Reynolds  Number  can  be  deter- 
mined from 


Re     = 

where  y     - 
h   = 


Q.  h/3 
V   RQZ 


(2.22) 


-  kinematic  viscosity 


thickness  of  vortex  chamber 


With  the  v^lue  of  Re  determined,  the  experimental  data  given 
by  Sarpkaya  on  viscous  efficiencies  can  be  entered  and  a  value 
of  E  determined. 

Knowing  all  the  values  of  the  variables  in  equation 
(2.20)  except  rD0,  allows  one  toi  solve  for  this  variable.   The 
position  of  the  probe  holes  from  the  center  of  the  sink-tube 
(distance  a  in  Figure  2.3)  can  then  be  determined  from  the 
relationship 
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r 
po 


--V^H^T  • 


which  can  be  verified  by  inspection  of  Figures  2.3  and  2.5. 
One  word  of  caution  must  be  made  concerning  a  minimum,  value 
of  rp0.   Viscosity  plays  an  increasingly  important  rol3  as 
r  approaches  zero  in  a  real  vortex.  If  rD0  is  too  large, 
full  advantage  is  not  made  of  the  amplifying  properties  of 
the  vortex.   On  the  other  hand,  if  it  is  made  too  small,  vis- 
cosity rapidly  decreases  the  strength  of  the  vortex.   This 
phenomena'  remains  to  be  investigated  more  fully  and  is  men- 
tioned here  only  for  design  considerations. 

The  design  data  used  in  constructing  the  rate  sensor 
is  summarized  in  Table  1.   The  material  used  in  construction 
was  plexiglass.   The  porous  ring  was  mace  of  felt  faced  with 
fine  wire  screening  which  gave  rigidity  to  the  ring  and  acted 
as  a  flow  straight ener.  The  probe  was  constructed  by  cutting 
a  1  3/4  inch  piece  of  number  19  gauge  stainless  steel  hypo- 
dermic tubing  in  half  and  inserting  a  plug.   The  number  80 
Drill  oize  holes  were  then  drilled  0.095  inches  apart.   This 
construction  allowed  the  holes  to  be  aligned  by  rotating  one 
side  of  the  probe  relative  to  the  other. 

2.3  Signal  amplification 

The  magnitude  of  the  signal  output  (pressure  and  flow) 
by  itself  is  too  small  to  be  used  conveniently  in  the  control 
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_;  TABLS  I 

w 

k  SIBJMARY  OF  RATE  SENSOH  DEoIGN  DATA 


Specified  Data: 

CO  max     Z       70°/sec 
Hq       =       3  in. 


A*max  =•  0.1     lbf/in' 


A# 


max 


Arbitrarily  Chosen  Data: 


•r3       -       0.1562  in. 


Derived  Data: 


0.0212  in.  (//19  gauge  hypodermic 
tubing  0D  x  .ID  s 
0.0425  x  0.0270  inj 


probe  hole  diameter  =  0.0135  in.  (^30  Drill) 

f 


r^  =  0.0512  in. 


a  -  0.0475  in. 

U3  =  171.2  ft/sec 

Q.  =  9.05  ft3/min 

S  =  0 .  50 
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system  loop  and  therefore  must  be  amplified.  The  state  of 
the  art  of  pure  fluic1  signal  amplification,  though  not  yet 
mature,  is  advanced  to  the  point  v;here  a  combination  of  ex- 
perimental and  analytical  methods  can  be  used  to  obtain  rea- 
sonable signal  amplification  with  tolerable  signal  to  noise 
ratios. 

The  design  problems  faced  at  this  stege  of  the  develop-  • 
'ment  of  the  rate  sensor  are: 

1)  Given  the  output  pressure-flow  relations  of  the 
rate  sensor  probe,  whet  are  the  supply  pressure  settings 
needed  at  each  amplifier  stege  to  insure  linearity  of  the 
signal? 

2)  How  many  stages  of  amplification  are  needed  to 
raise  the  pressure  and  f low> to  acceptable  levels? 

The  last  question  could  not  be  answered  at  this 
stage  of  the  design  for  the  rest  of  the  systems  com- 
ponents had  not  yet  been  designed.   Therefore,  it  was 
decided  to  use  two  stages  of  amplification  and  after 
the  amplifier  stages  were  attached,  to  accept  the  out- 
put pressure-flow  relations  of  the  final  amplification 
stege  as  the  output  characteristics  of  the  sensor. 

The  first  question,  although  a  relatively  straight  for- 
ward design  Question,  is  covered  only  lightly  in  the  literature 
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The  procedure  used  to  obtain  the  supply  pressure  settings 
described  belcw  is  by  no  means  exact  but  it  will  let  the 
designer  set  the  pressures  close  to  their  optimum  settings 
on  the  first  try. 

The  amplifiers  used  to  amplify  the  rate  signal  output 
were  Corning  Glass  Work's -Standard  Size  Proportional  Amplifiers. 
The  important  dimensions  of  these  devices  are: 

power  nozzle  width  ■:  w  s.  0.010  in. 
control  nozzle  width  s  1.5w 

depth  -  2.5w 
•  receiver  distance  -  lOw 
Figure  2.7  is  a  sketch  of  the  interaction  region  of  a  typi- 
cal proportional  amplifier  which  operates  principly  by  mo- 
mentum interaction.  Following  the  development  given  by 

5 
Dexter  ,  the  deflection  of  the  power  stream  due  to  one  control 

momentum  only  is  determined  as  follows: 

tan  Q     - 


H 


Va 


where   V.        r  mc   Vc 


mo+mp 

\   ■-    hs 

mc+mp 
and   mp  *  mass  flow  rate  of  the  power  jet 

« 

# 

cic  -  mess  flow  rate  of  the  control  jet 
\ 

{  Vc  =  velocity  of  the  control'  jet 

Vp  s velocity  of  the  power  jet 
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-  nozzle 


Lgure  2.7     ' 
Typical  Proportion-1!   Pure  Fluid  Anvolifier 


•h 
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\   tan  Q   =   mc  vc 

*P  vp 

but  mv  =  p^V2 

,\  tan   -  P*c  V2   • 

:   Pdwc  v02 

PA^  V  2 
r  P  p 

P  d  wp  Vp2 

where   d  r  depth 

wc  vc2 

2//0 

1 

wp  VP2 

p 

wc  2/P  Pc 

=  wo  po 

WP  ~'r  *p 

WP  pp 

wc  =  control  nozzle  width 
Wp  r  power  nozzle  width 

.'.   tan#  =         "e  ^  (3.24) 

but  V2  = 
/.  tan£  =    w°  "'f fc   =  w0  P„      12-.2SJ 

wp  27^p      ^rv 

He illy  and  Moynihan18  indicate  that  the  output  of  a 

proportional  device  will  remain  linear  provided  the  deflection 

of  the  power  jet  at  the  receiver  ports  does  not  exceed  one 

half  the  nozzle  width.  Thus: 

tan  Q^- wp/2      for  linearity 

id 

where  L  =  receiver  distance 

Combining  this  with  equation  (2.25)  gives: 

w-   ^  wc  pc    for  linearity 


2L      wp  Pp 
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This  result  is  readily  extended  to  two  control  ports  with  the 
result: 

WP  *    wc  (^cL-Pch)    for        (2'26) 


2L        v/p  P       linearity 

Having  developed  the  basic  equations  of  the  amplifier, 
the  output  characteristics  of  the  rate  sensor  probe  must  now 
be  investigated.  According  to  potential  flow  theory,  if  the 
probe  is  set  so  that  the  probe  holes  are  at  45  degree  angles  to 
the  oncoming  flow,  the  pressure  in  the  probe  chambers  is.  below 
ambient.  Since  the  probe  chambers  are  connected  to  the  control 
ports  of  the  first  stage  amplifier,  this  means  that  only  the 
pressure  differential  existing  in  the  two  chambers  and  com- 
municated to  the  first  stage  control  ports  is  available  for 
deflecting  the  jet.  This  is  undesirable  as  it  results  in 
less  gain  when  compared  to  operation  with  positive  flow  through 
both  control  ports.  However,  due  to  viscosity  the  ideal  flow 
pattern  is  radically  altered  so  that  the  pressure  existing  in 
the  chambers  when  the  flow  angle  is  set  to  45  degrees  is  posi- 
tive. 

The  various  flow  rates  to  be  expected  out  of  either 
end  of  the  probe,  as  the  main  sink  flow  is  deflected,  can  be 
approximated  by  setting  the  sink  flow  at  its  operating  value 
and  turning  the  probe  with  respect  to  the  flow.  The  results 
of  such  an  experiment  are  shown  in  figure  2.8. 

The  output  pressure-flow  characteristics  for  a  given 
?b  are  shown  in  Figure  2.9.  An  infinite  number  of  such  curves 
exist  since  Pb  varies  with  A  Q* 
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Pb  (in.  alcohol) 
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Figure  2.3 
Probe  Boundary  Pressure  versus  Mo  T.oad  Flow 
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Figure  2.9 
Probe  Outnut  Pressure*- Flow  Characteristics 
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i^      .  p  ^)      ■  =  A p 

ucL   ^cii'max    ^  c 


The  procedure  initially  used  to  set  the  first  stage 
supply  pressure  (Py  )  was  to  determine  the  output  pressure-  ' 

pJL 

flow  relationship  for  maximum  A  P  and  superimpose  the  input 
characteristics  of  the  first  stage  amplifier  to  determine 

max*  ?hofl>  from  equation  (2.25),  the 
supply  pressure  was  determined.  However,  when  this  procedure 
was  followed,  it  was  found  that,  because  the  control  port  area 
of  the  first  stage  amplifier  was  so  large  compared  to  the  probe 
hole  area  (nearly  three  times  as  large),  the  input  impedance 
of  the  amplifier  was  nearly  zero. 

Circumventing  this  problem,  equation  (2.24)  was  used 
in  lieu  of  (2.25).  With  the  constraint  placed  on  tan  Q  and 
for  two  control  port  operation,  equation  (2.24)  becomes: 

w?     .  "c  <ycl2  -  vca8W  <a.a?) 

2L        wp  Vp2 

From  measurements  of  the  flow  rate  occurring  at  the 

expected  maximum  pressure  just  outside  the  probe  hole, 

(V  T2  -  V  r>2)     was  obtained  from  V~L2  -  V_^2  ■  1 
v  cL     cR  max  CJ-     0£i  ,  £ 

^cL2max  "  ^cH2min^  wilich  then  allowed  Vp  to  be  obtained  from 

equation  (2.27) 

To  obtain  Pp  the  following  relationship  was  used: 

%  ■-  cdP    ap    V*/P  ?P  (2-28) 

where  Qp  s  Ap  Vp 

C,  =  discharge  coeficient  of  the  power 
^   nozzle  which  was  determined  from  the 
design  data  sheet  published  by  Corning 
glass  '.,'orks. 
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From  the  calculations  in  Appendix  A. 2,  r  .  was  determined  to 
be  7.7  inches  of  water. 

The  second  stage  pressure  was  set  using  information 
given  in  the  design  data  sheet  published  by  Corning  Glass 
Works  for  these  amplifiers.  This  sheet  states  that  for  max- 
imum gain  the  control  bias  level  pressure  shoulc  be  one  tenth 
the  power  jet  pressure*. 

The  bias  level  of  the  second  stage  controls  was  deter- 
mined by  setting  the  first  stage  supply  pressure  at  7.7 
inches  of  water  and  measuring  the  pressure  in  the  line  connect- 
ing the  output  of  the  first  stage  and  the  control  port  of  the 
first  stage.  This  procedure  resulted  in  an  initial  setting  of 
Pp2  equal  to  5  psi. 

2.4  Hate  Sensor  Dynamics 

This  rate  sensor  is  characterized  dynamically  by  1) 
a  pure  time  delay  resulting  from  the  time  it  takes  for  the  flu- 
id particles  containing  tangential  momentum  to  reach  the  probe, 
2)  a  first  order  lag  due  to  the  resistances  of  the  probe  holes 
and  control  nozzles  of  the  first  stage  amplifier  and  the  vol- 
ume connecting  them,  and  3)  the  dynamics  of  the  staging  am-  • 
plifiers  and  the  lines  used  in  interconnection. 

The  transport  delay  is  the  most  significant.   It  may 
i 


This  procedure  used  in  setting  the  first  stage 
resu 

was  not  used. 


pressure  resulted  in  a  P  ,  ■  1.95  inches  of  water  and  hence 
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be  approximated  by  the  tine  it  takes  for  the  sensor  to  com- 
pletely  renew  the  supply  of  fluid  in  the  vortex  chamber 
and  sink-tubes.  Thus: 


7"  =   — - —   seconds  (2.29) 

where  V  -  total  volume 

Q,  =  flow  rate  given  by  equation  {2.23) 

The  total  volume  (V)  is  given  by: 

V  =   <TT-H02h  +  2^rrslp 

where  H0 ,  h  and  rs  are  as  before 


1  z     distance  to  probe  measured  along  the 
P    axis  of  the  sink-tube 


For  the  geometry  and  flow  rates  of  this  device,  the  value  of 
T  is  0.0556  seconds. 

The  time  constant  of  the  first  order  lag  referred  to 
in  item  2)  was  found  to  be  extremely  small  and  was  neglected 
in  arriving  at  the  transfer  function  of  the  rate  sensor. 

The  dynamics  of  the  staging  amplifiers  are  discussed  by 
Belsterling  and  Tsuis  and  Boothe7.  These  works  show  the  de- 
vices to  respond  very  much  faster  than  the  time  delay  of  this 
rate  sensor.   In  addition,  the  design  data  sheet  for  these 
amplifiers  indicate  that  the  frequency  gain  is  over  1000  cps, 
well  out  of  the  range' of  the  rate  sensor. 

Thus  the  transfer  function  relating  the  output  (A?) 
to  the  input  angular  rate  (CJJ  is  given  by: 
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Cu  x 

where  •  K^    =  steady  stete  gain 

» 

Tt    -  pure  tine  delay  given  by 
equation  ( 29) 

i 

S     -  Laplace  operator  d/dt 

2.5  Testing  Procedure  and  Results 

The  rate  sensor  was  made  into  a  self § contained  pack- 
age by  attaching  the  staging  amplifiers  directly  to  one  side 
of  the  vortex  chamber.  The  sensor  was  equipped  with  a  man- 
ifold allowing  air  to  be  supplied  with  only  one  line.  This 
manifold  pressure  was  then  reduced  to  the  required  pressures 
needed  at  the  vortex  chamber  supply  plenum  and  power  jets 
of  the  staging  amplifiers.   Variable  resistors  were  used  to 
set  the  pressures  and  were  constructed  from  tapered  dowel 
pins  through  which  holes  were  drilled.   Photographs  of  the 
assembled  rate  sensor  are  shown  in  figures  7.4 and  D.l  .   A 
differential  pressure  to  electric  transducer  was  attached  to 
the  outputs  of  the  second  stage  amplifier  and  the  assembled 
device  was  tested  on  a  rate  table  at  the- Instrumentation 
Laboratory  of  the  Massachusetts  Institute  of  Technology. 
The  testing  procedure  consisted  of  the  following: 
1')   Steady  state  test  over  the  range  1.10  £  60^-1.15 
radians  per  second  with  the  output  loaded  by  orifices  approx- 
imating the  size  of  the  control  parts  of  an  amplifier  similiar 
to  those  used  on  the  rate  sensor  for  signal  amplification.  • 


3>» 


i 
2}   Steady  state -test  over  the  same  range  with  the  out- 

put  flow  equal  to  zero;  i.e.,  infinite  load. 

3)   frequency  response  test  at  discrete  frequencies 
up  to  4.67  cps. 

4}   Steady  state  test  while  the  device  was  being  simuV 
taneously  rotated  about  its  sensing  ax^s  and  an  arbitrary  axis. 

The  test  results  for  the  first  three  tests  are  shown  in 
Figures  2.10,  2.11,  2.12,  and  2.13. 

2.6  Discussion  of  Results 

The  results  of  the  steady  state  rate  tests  show  the 
following: 

1.  The  very  small  signals  available  from  the  probe  can 
indeed  be  amplified  into  reasonable  signals. 

2.  The  output  of  the  sensor  is  not  quite  linear  over  — 
the  range  tested. 

3.  The  gain  of  the  sensor  with  clockwise  rotation  is 
significantly  larger  than  the  gain  with  counter  clockwise  ro-* 
tat  ion  when  the  bias  is  subtracted  out  of  the  results. 

4.  The  gain  increases  continually  as  the  input  rate 
is  increased. 

5.  The  dog  leg  effect  produces  a  highly  undesirable 
rapid  change  of  gain  in  the  counter  clockwise  direction. 

In  an  attempt  to 'explain  the  last  four  items,  accurate 
flow  ajid  pressure  measurements  were  taken  on  each  side  of 
the  probe  as  the  angle  of  incidence  of  the  probe  holes  with 
the  oncoming  flow  was  vsried.  Originally,  in  aligning  the 
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Figure  2.13 
Typical  Vortex  Rate  Sensor  Frequency  Response  Data 


39 


I 

probe  holes  (the  turning  of  one  probe  hole  with  the  other  held 
fixed),  only  the  dead-ended  pressures  were  used,  talcing  the 
probe  holes  as  being  aligned  when  the  pressures  were  equal. 

As  was  pointed  out  in  article  2,3,  signal  amplifica- 
tion is  essentially  dependent  on  the  flows  cooing  from  the 
probes.  When  these  flows  were  measured,  it  was  found  that 
equal  dead-ended  probe  pressures  resulted  in  radically  dif- 
ferent flows  indicating  that,  although  the  holes  were  sup- 
posedly of  the  same  area,  they  were  behaving  as  if  they  were 
not. 

The  probe  was  removed  from  the  sensor  and  observed  under 
a  50  power  microsoope.  This  showed  that  during  the  machining 
process,  the  drill  had  been  allowed  to  out  a  trough  into  the 
entrance  to  one  of  the  holes.  The  effect  of  this  trough  was 
apparently  to  catch  more  flow  from  the  main  stream  resulting 
in  different  orifice  characteristics. 

The  effect  of  having  the  two  probe  holes  producing 
different  flows  for  equal  boundary  pressures  was  as  follows: 
It  has  previously  been  shown  that,  as  the  sink-tube  axial 
velocity  is  deflected,  the  pressure  outside  the  probe  holes 
increases  or  decreases  depending  on  the  sense  of  angular  ro- 
tation and  the  probe  hole  in  question.  As  the  pressure  out- 
side the  probe  hole  which  gives  less  flow  for  a  given  pressure 
is  reduced,  a  point  is  reached  where  the  flow  approaches  zero. 
By  referring  to  Figure  2.8,  it  can  be  seen  that,  as  the  flow 
approaches  zero,  the  relation  between  the  pressure  outside 
the  hole  and  the  resulting  flow  beoomes  highly  non-linear. 
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This  highly  non-linear  flow  region  produces  the  dogleg  ef- 
fect noted  above.   This  effect  is  not  present  when  the  direc- 
tion of  rotation  is  reversed  for,  as  the  boundary  pressure 
outside  the  hole  which  gives  the  greater  flow  for  a  given 
pressure  is  reduced,  the  non-linear  portion  of  the  curve  is 
avoided*.  At  the  same  time,  the  other  probe  hole  is  experi- 
encing increasing  pressure  and  hence  this  combination  produces 
no  dogleg  effect.  Thus,  to  remove  the  effect,  the  probe  hole 
without  the  trough  was  altered  to  give  approximately  the  same 
pressure-flow  characteristics  as  the  other  probe  hole. 

To  explain  the  difference  in  the  maximum  signal  output 
between  clockwise  and  counter  clockwise  rotation,  one  has  only 
to  realize  that  the  positioning  of  the  probe  holes  with  respect 
to  the  oenter  line  of  the  sink-tube  is  critical.   If  the  holes 
are  off  center  the  slightest  amount,  the  output  will  be  de- 
creased on  one  side  and  increased  on  the  other  resulting  in 
different  peak  outputs. 

The  effect  of  the  gain  continually  increasing  with  in- 
creasing CJcan  be  partly  explained  by  realizing  that  deflec- 
tion of  the  first  stage  power  jet  is  roughly  proportional  to 
(A«i)2.  'Generally  the  output  pressure  of  a  proportional  am- 
plifier is  given  by 

Ap0  =  kxAp0 

but  Ap0  =   K2  (A*2) 


At  least  over  the  range  of  rotations  tested. 
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hence  A?Q   -  K5  (A<J2 

Now  as  shown  in  Figure  2.8,  A**  is  approximately  proportional 
to  CO  over  a  substantial  range  of  operation,  making  the  output 
pressure  generally  proportional  to  the  square  of  6J  .   This 
of  course  is  only  approximate  but  does  roughly  explain  the 
increasing  gain  of.  the  device.   In  fact,  if  the  output  pres- 
sure  is  converted  to  flow  using  the  output  characteristics 
of  the  second  stage  amplifier,  one  finds  the  output  flow  quite 
linear  with  respect  to  U) • 

Although  accurate  verifications  of  the  explanations 
given  for  items  2,  3,  4,  and  5  above  were  not  accomplished 
due  to  inaccessibility  of  the  rate  table,  the  explanations 
g^ven  were  partially  verified  and  showed  that  the  dogleg  had 
been  removed  but  that  only  slight  improvement  was  made  in 
balancing  the  outputs  by  centering  the  probe. 

The  zero  input  bias  signal  was  found  to  be  only  partly 
due  to  the  unequal  flows  from  the  probe.  The  amplifiers 
themselves  are  biased  when  the  control  flows  are  zero.  This 
can  be  reduced  somewhat  by  proper  arrangement  of  the  cascade 
so  that  the  output  bias  of  one  tends  to  cancel  the  output 
bias  of  the  next  stage.  However,  adjusting  orifices  general- 
ly must  be  placed  somewhere  in  the  cascade  to  remove  the  bias. 

« 

The  results  of  the  frequency  response  tests  verified  the 
pure  time  delay  nature  of  the  dynamics  of  this  sensor.   The 
experimental  phase  data  is  approximate  at  the  low  frequencies 
because  the  method  used  to  obtain  the  data  (the  input  and 
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output  were  recorded  on  permanent  recording  paper  and  were 
compared  in  order-  to  measure  phasej  could  not  discern  small 
phase  shifts.  The  amplitude  data,  however,  is  quite  accurate, 

Lissajous  figures  on  an  osciloscope  could  not  be  used 
to  obtain  frequency  response*  data  because  the  output  noise 
moved  the  figures  around  the  scope  precluding  any  measurements*. 

The  best  data  was  obtained  by  taking  pictures  of  the 
oscilloscope  traces  and  measuring  phase  and  attentuation  data 

from  these.  The  rapid  increase  in  phase  verifies  the  existence 
of  a  pure  time  delay.  The  deviation  from  theory  in  phase 
shift  and  attenuation  is  probably  due  to  volumes  in  the  lines 
connecting  the  first  and  second  stage  amplifiers  and  the  vol- 
ume and  resistance  of  the  orifices  between  the  output  of  the 
second  stage  amplifier  and  the  pressure  transducer.      , 

The  results  of  the  last  test  showed  the  output  signal 

did  not  change  while  the  sensor  was  being  simultaneously 

i 
rotated  about  its  sensing  axis  and  an  arbitrary  axis. 

2.7  Conclusions  and  Recommendations 

Throughout  the  development  of  this  rate  sensor  little 
attempt  has  been  made  at  optimization.   The  objective  was  to 
develop  a  rate  sensor  capable  of  being  used  in  a  control  system. 
To  this  end  the  following  conclusions  are  drawn: 


*  The  noise  of  fluicic  proportional  amplifiers  is  of 
a  peculiar  nature.  Although  high  frequency  noise  is  prevalent 
as  shown  in  Figure  2.13,  a  more  disturbing  very  low  frequency 
noise  is  present  which  causes  the  signal  to  wander  and  make 
unexpected  excursions  which  appear  almost  as  steady  state  bias 
changes. 
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1.  Use  of  the  design  procedure  suggested  in  article  2.2 
will  allow  the  designer  to  theoretically  design  a  rate  sensor 
peculiar  to  his  applications. 

2.  The  steady  state  gain  of  the  unamplified  signal 
may  be  predicted  analytically  but  due  to  the  low  input  im- 
pedance of  the  type  of  amplifier  used,  the  amplified  gain  will 
have  to  be  determined  experimentally. 

3.  The  dynamics  of  the  sensor  may  be  determined  to  a 
reasonable  degree  of  accuracy  by  theory. 

4.  The  output  of  the  sensor  is  very  dependent  on  the 
accuracy  of  the  machining  of  the  probe  hole3  and  alignment  in 
the  sink-tube. 

5.  Noise  of  the  very  low  frequency  type  is  a  serious 
problem  causing  unexpected  wandering  of  the  DC  level  of  the 
signals. 

6.  Contrary  to  the  results  of  Sarpkaya,  it  was  found 
that  closing  one  of  the  sink-tubes  and  keeping  the  flow  con- 
stant, resulted  in  an  increased  output  signal  for  a  given  CO  • 
•Since  this  equivalently  means  nearly  the  same  signal  with  less 
power  consumption,  the  author  recommends  that  sensors  of  this 
type  be  built  with  only  one  sink-tube  and,  as  recommended  by 
Sarpkaya,  that  the  end  opposite  the  sink-tube  be  designed  to 
streamline  the  flow  into  the  sink.-tube. 
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CHAPTER  III 
ANGULAR  POSITION  SENSOR  ANALYSIS  AND  DESIGN 

3.1  Criteria  Used  in  Selecting  Operating,  Principle 

Angular  position  of  the  system  can  be  obtained  by 
analog  or  digital3'9  integration  of  the  angular  rate  signal 
available  from  the  vortex  rate  sensor.   Digital  integration 
would  possibly  involve  developing  a  pressure  sensitive  os- 
cillator and  some  form  of  digital-to-analog  converter,  both 
of  which  are  beyond  the  time  limitations  of  this  thesis. 
Since  a  linear  position  signal  over  a  large  range  of  angular 
positions  for  demonstration  purposes  was  desired,  analog 
integration  was  ruled  out  since  very  large  'volumes  would  have 

* 

been  required  toustore"the  position  sig£al.   furthermore,  the 

7 

state  of  the  art  of  pure  fluid  analog  computation  is  such 
that  up-down  integration  is  inaccurate.        ; 

3.2  Position  Sensor  Description 

Because  of  the  problems  mentioned  above,  it  was  decided 

•i 

to  use  a  moving  part  position  sensor  consisting  of  a  pendulum 
set  between  a  p8ir  of  flapper  nozzle  valves  to  obtain  a  dif- 
ferential pressure  signal  that  is  a  function  of  angular  posi- 
tion. A  schematic  of  this  sensor  is  shown  in  Figure  3.1  and 
a  photograph  of  the  final  design  is  shown  in  Figure  D.3.   The 
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device  is  designed  so  that  as  the  sensor  rotates,  the  pendulum 
trys  to  return  to  a  vertical  position.   In  so  doing,  it  widens 
the  gap  on  one  side  anc  closes  the  gap  on  the  other  until  the 
moments  caused  by  the  flow  forces  from  the  nozzles,  balance 
the  moment  due  to  the  weight  of  the  pendulum.   .Vhen  this 
occurs  a  differential  pressure  exists  in  the  chambers  and 
constitutes  the  output  signal  of  the  sensor. 

It  is  readily  seen  that  the  output  differential  pres- 
sure signal  varies  as  the  sine  of  the  angular  position  and 
that  the  sensor  also  senses  direction  of  angular  position. 

3.3  Discussion  of  the  Dynamics  of  a  Pendulum  Angular 
Position  Sensor 

If  the  bearing  which  supports  the  pendulum  is  relative- 
ly free  of  damping,  the  angular  position  of  the  pendulum  is 
unstable  and  the  slightest  perturbation  will  start  the  device 
oscillating.   In  the  bang-oang  control  envisioned  for  the 
final  system,  this  oscillation  is  undesirable.   To  understand 

the  reasons  for  this  instability,  a  linear  model  of  the  system 

f 
can  be  constructed  to  describe  the  behavior  of  the  system 

when  the  angular  position  of  the  pendulum  is  small.   Figure  3.2 

shows  an  analytical  model  of  the  sensor.  .Summing  moments 

about  the  pivot  point  gives  the  following  differential  equa-  • 

tion: 
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where  Q        =     d^/dt 

bl^  =  equivalent  viscous  damping 


k    -   equivalent  spring  constant  of  the 
flapper  nozzle  valves 


Tj   =  a  disturbance  torque 

Neglecting  the  term  \lQ   ,  since  it  is  much  smaller  than 

W 


[ 


Js2  +  B.s  +  X 


kln#,  and  introducing  the  symbols  J*  -JS-  lw  ,   B=  bl^, 
K=  kln,  and  s=-  d/dt  gives  • 

- 

9     s  ?d         (3-D 

The  spring  force  K  $is  given  by 

KQ     =   (PiA  -  P2A)ln    =   APAla 
where  P^_   =  pressure  in  one  chamber 

?2   ■  pressure  in  the  opposite  chamber 
rearranging  gives 

X    =  Aln  A?  (3.2) 

9 

The  linearized  transfer  function  between  AP  and$ 
is,  in  general,  a  first  order  lag  of  the  form  . 

AP    -    C1  (3.3) 

9  "Ts+1 
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Substituting  equation  (3.3)  into  (3,2)  gives 


AlnC 


n°l 


C2 


(3.4)- 


Tsvl      Ts+1 
If  this  equation  is  substituted  into  equation  (3.1)  the 
result  is 


[ 


Js2  +  Bs  +  £2. 


B   =   * 


Ts+1  J 
From  this  the  transfer  function  between  Q  and  T^  is 


j£ 


Ts  +  1 


JTs^(J+BT)s2f  bs>C2 


The  characteristic  equation  of  this  transfer  function  is 


TjTs3  +  (J  +  BT)  s2  4-  £3  ♦  C21 


The  Routh  stability  criterion  for  linear  systems  gives  the 
condition  for  stability  as 

B(J>BT)      >       ZTQo 


or 


B^VJ2  ♦  4JT^C2         -J 


(3.5) 
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where  the  positive  sign  for  the  radical  has  been  chosen  since 
the  damping  can  not  be  negative.  If  equation  (3.5)  is  re- 
written as 


B 


•  «[Vj?p 


4JC2 


jr. 

T 
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it  can  be  seen  that  the  value  of  C2  is  the  important  para- 
meter determining  the  stability  of  the  system.  As  Cg  in- 
creases, 3  must  be  made  correspondingly  larger.  ]?rom 
equation  (3.4),  C2  is  the  value  of  the  steady  state  spring 
constant  of  the  flapper-nozzle  valves.   Thus,  as  this  spring 
constant  increases,  the  magnitude  of  the  damping  must  be 
increased. 

Initially,  to  introduce  damping  into  the  system,  the 
ball  bearing  used  to  support  the  pendulum  was  packed  with 
heavy  grease.  This  aid  not  introduce  a  sufficient  amount 
of  damping  into  the  system.  The  problem  was  finally  solved 
by  placing  two  small  squeeze  film  dampers  on  either  side  of 
the  pendulum  snd  adding  a  minute  amount  of  grease  to  aid  in 
the  damping  process  . 

3.4  The  Influence  of  Acceleration  on  the  Output  Signal 

The  output  signal  of  this  sensor  will  be  influenced 
by  the  angular  acceleration  of  the  device  to  which  it  is 
attached.  Figure  3.3  is  a  sketch  showing  the  position  sen- 
sor mounted  on  the  system  and  offset  from  the  axis  of  rotation 
by  a  distance  r. 

If  the  sensor  is  given  an  angular  acceleration  &  , 


*  Over  most  pressure  settings  the  air  in  the  squeeze 
film  dampers  provided  sufficient  damping.  However,  if  the 
chamber  pressures  were  made  too  high  then  the  size  dampers 
used  were  inadequate. 
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then  the  end  of  the  pendulum  receives  a  linear  acceleration 
rd.  The  torque  on  the  pendulum  is  given  by 

T^   -    w     ra  (3.6) 

g 

.  This  torque  must  be  balanced  by  the  pressure  forces  on  the 
end  of  the  pendulum  giving 

Apli  Aln  =   .».  ra 
or 

This  A?  corresponds  to  an  error  signal.  It  may  be  reduced 
by  making  r  small  assuming  the  other  parameters  are  fixed. 

The  moment  of  inertia  (1)  of  the  pendulum  is  given  by 

I   =  -2L  1  2 

g   w 

Thus  the  resulting  torque  on  the  pendulum  is  given  by 

I  TAA   =  If  Xw2  0.  (3.8) 

This  torque  too  must  be  balanced  by  the  pressure  forces  on  the 
end  of  the  pendulum  giving 


APAA  *1n 

w 

g 

1  2 
*w   a 

ApAA     ■  = 

Wu.    w 

a 

or 

.....  2 

(3.9) 
6Aln 


*  Subscripts  LA  end  AA  refer  to  linear  acceleration 
anc"  angular  acceleration  respectively. 
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This  A  F  also  represents  an  error  in  measuring  angular  position. 
The  total  error  signal  is  the  sum  of  the  two  differential  pres- 
sures given  by  equations  (3.7)  and  3.9) 

A?«**'^rif7-   f-fe-)a  (3>10) 

This  pressure  error  may  be  converted  into  an  equivalent  angu- 
lar position  error  by  the  equation 

# 

d   =  Sin"1  (A?  Aln  )  (3.n) 

\   wlw    / 
I 

which  is  obtained  by  a  steady  state  torque  balance  on  the 

•i  _ 

pendulum  when  at  an  angle^  from  the  vertical.   Substitution 

of  the  value  of  A?  from  equation  (3.10)  into  equation  (3.11) 

will  give  the  error  in  terms  of  angular  position.   Primarily 

because  the  accelerations  to  be  cea&t  with  in  this  system  are 

quite  small,  substitution  of  the  actual  values  into  equations 

(3.10)  and  (3.11)  showed  that  the  maximum  probable  error  v/ould 

be  less  than  one  degree.  If  this  system  were  to  be  made  into 

a  very  fast  high  torque  system,  this  type  of  position  sensor 

woulfi  probably  be  undesirable  unless  the  designer  wished  to 

have  acceleration  feed  back  in  the  system. 

3.5  Transfer  Function 

Because  the  amount  of  damping  added  to  the  system  to 
make  it  stable  is  difficult. to  determine,  it  was  deoided  to 
connect  the  output  signal  into  a  pressure  to  electric  trans- 
ducer and  experimentally  determine  the  transfer  function.   By 
observing  step  responses  on  an  oscilloscope,  the  system  dynamics 
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can  readily  be  determined. 

The  procedure  was  as  follows: 

V 

1)  The  output  signal  was  connected  to  the  transducer 
and  loaded  into  orifices  which  approximated  the  size  to  be 
expected  at  the  control  ports  of  the  summing  device. 

2)  The  output  signal  was  then  displayed  on  an  oscil- 
loscope. 

3)  The  pendulum  was  given  an  initial  deflection  and 
then  released. 

The  response  of  the  pendulum  was  observed  to  be  es- 
sentially that  of  a  first  order  lag.   3y  measuring  the  time  it 
took  the  signal  to  decay  to  36.8  percent  of  Its  initial  value, 
gave  a  close  approximation  to  the  time  constant  of  the  sensor. 
By  averaging  a  large  number  of  such  measurements  the  time  con- 
stant of  the  system  was  found  to  be  0.04  seconds.  This  time 
constant  was  found  to  be  much  larger  than  the  time  constant 
associated  with  the  flapper-nozzle  chambers*.  Apparently, 
the  measured  time  constant  is  essentially  that  of  the  mass- 
damping  time  constant  of  the  system. 

A  typical  steady  state  output  curve  is  shown  in 
Figure  3.4  from  which  the  steady  state  gain  may  be  measured. 

The- transfer  function,  relating  the  output  pressure 
signal  to  the  input  angular  position  signal  for  small  valves, 
of  Q  ,  is  given  by 


*  The  time  constant  cue  to  the  volume  of  the  flapper 
nozzle  chambers  and  connecting  lines  was  found  to  be  0.008 
seconds.  The  calculations  are  not  shown  here  as  they  can  be 
found- in  any  standard  reference  such  as  reference  10. 
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A?    =       K2  (3.12) 

0  7V>fl 

. 

where  K2   -     gain  to  be  determined  from  a  plot 

similar  to  figure  3.4  once  the  final 
system  pressures  have  been  set 

7"2   ~  0.04  seconds 

Before  concluding  this  chapter  it  is  to  be  noticed 
that  little  information  has  been  given  concerning  design 
procedures.  This  was  purposely  done  as  such  procedures  are 
more  adequately  covered  in  the  literature^-0.  However,  the 
designer  must  keep  in  mind  the  over  all  system  and  the  re- 
lationship of  each  component  to  the  others.   For  example, 
it  is  obvious  that  the  output  of  the  position  signal  is  to 
eventually  be  compared  with  the  rate  signal.  It  is  therefore 
reasonable  to  assume  that  these  signals  should  be  of  the 
same  order  of  magnitude. 

3.6  Switching  Curve  Computer 

As  shown  in  Appendix  0.1,  the  ideal  switching  curve 
for  this  system  is  a  parabola  passing  through  the  origin. 
The  slope  of  this  curve  is  infinite  at  the  origin.   Since  it 
was  planned  to  generate  the  switching  curve  by  operating  on 
the  output  signal  of  the  position  sensor,  this  meant  that  as 
Q  varied  slightly  from  its  zero  position  an  infinite  gain 
would  be  required.   Because  of  this  consideration  and  the 
fact  that  the  output  signal  of  the  rate  sensor  is  not  depend- 
able at  very  low  angular  rates,  the  output  of  the  position 
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sensor,  without  modification,  was  used  as  the  switching  curve. 
Such  a  sine  curve  output  can  be  adjusted  to  approximate 
a  parabola  over  a  large  range  of  position  signals  as  shown  in 
Figure  3,4,  The  approximation  is  poor  when  Q  is  small,  how- 
ever, and  thus,  it  is  to  be  expected  that  "linear  switching" 
will  take  place  when  the  errors  introduced  into  the  system 
are  small. 


v 
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CHAPTER  IV 
C 

SIGNAL  COMPARATOR 


4.1  General 

The  method   selected  to  algebraically  sum  the  pressure 
signals  from  the  rate  sensor  and  position  sensor  consists  of 
adding  additional  sets  of  controls  to  a  proportional  amplifier 
and  having  the  power  jet  deflect  in  accordance  with  the  vector 
sum  of  their  momentums*.  A  stream  interaction  proportional 
-amplifier,  operating  on  this  principle,  is  available  from 
Corning  Glass  Works**  and  has  two  sets  of  control  pox'ts.   The 
design  data  sheet  accompanying  this  amplifier  does  not  in- 
dicate that  *the  second  set  of  controls  can  be  used  for  sum- 
ming and,  therefore,  gives  no  data  on  the  gain  of  these  ports. 
In  fact,  the  ports  are  not  equipped  with  fittings  and,  there- 
fore, had  to  be  adapted  for  use  in  this  oontrol  system.  A 
photograph  of  this  "sunimer-amplif  ier"  is  shown  in  Figure  D  .2. 
The  vertical  controls  are  the  seconoary  controls  and  were  ■' 
adapted  for  use  by  inserting  hypodermic  tubing  into  the  holes 
and  then  cementing  the  tubing  into  place. 


*  Other  methods  of  summing  signals,  using  no  moving 
parts,  are  available  and  could  have  been  used  in  this  system. 

**  Catalog  Number  FD2511-2-1211  September,  1965 
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4.2  Experimental  Analysis 

A  silhouett.e  of  the  amplifier  described  above  is  given 
in  the  design  data  sheet.  From  this  silhouette,  it  was  found 
that  the  width  of  the  secondary  control  ports  is  less  than 
the  width  of  the  primary  ports.   In  eddition,  the  secondary 
ports  are  not  perpendicular  to  the  power  jet.   Because  of  these 
geometrical  differences,  it  is  obvious  that  if  the  same  pres- 
sure differential  is  applied  to  both  sets  of  controls  ports, 
the  power  jet  will  be  deflected  a  greater  amount  by  the  pri- 
mary control  ports.  As  a  consequence  of  this  greater  deflec- 
tion, the  output  pressure  will  be  larger.   Therefore,  the 
primary  control  ports  will  exhibit  a  greater  gain. 

To  determine,  the  difference  in  the  gain  of  these  two 
sets  of  control  nozzles,  a  number  of  tests  were  performed  on 
the  amplifier.   The  tests  consisted  of  1)  systematic  varia- 
tion of  the  control  port  bias  levels  about  the  levels  ex- 
' pected  from  the  outputs  of  the  rate  and  position  sensors, 
2)  systematic  variation  of  the  supply  pressure,  and  3)  at  any 
particular  bias  level' and  supply  pressure  setting,  recording 
the  gains  of  the  two  sets  of  control  ports. 

Typical  data  from  these  tests'  for  one  particular  com- 
bination of  supply  pressure  and  bias  level  settings  is  shown 
in  Figure  4.1  During  all  these  tests,  the  output  of  the 
summer-amplifier  was  loaded  into  a  manometer. 

4.3  Discussion  of  Results  and  Conclusions 

The  results  of  these  tests  showed  that  the  gain  ratio 
of  the  two  pairs  of  control  ports  did  not  remain  constant  as 
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Figure  l .  1 
^ical  Summer- Amplifier  Test  Results 
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the  control  bias  level  or  as  the  supply  pressure  was  changed. 
However,  for  any  given  set  of  control  bias  pressures  and  a 
given  supply  pressure,  the  ratio  of  the  gains  of  the  two  con- 
trol pox-ts  could  repeatedly  be  obtained  after  signals  were 
introduced  and  then  removed. 

The  results  further  showed  tfhat^  once  the  ratio  of  gains 
had  been  determined  at  a  particulars  setting  of  bias  levels  and 
supply  pressure  that  the  algebraic  sum  of  two  control  signals 
could  accurately  be  predicted. 

Prom  the  results  of  these  tests,  it  was  conoluded  that 
this  device  could  be  used  as  a  summing  device  with  accurate 
and  predictable  results  if  the  bias  levels  dt   the  input  sig- 
nals and  the  power  'jet  pressure  were  known,  - 
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CHAPTER  V 


i 


DRIVING;  BISTABLE  AMPLIFIED  WITH  PliOPOitTIONAL  AivIPLIPlERS 


5.1  Introduction  to, Input  Control  Impedances  of  Bistable 
Amplif iersx 

In  Figure  5.1  a  typical  bistable  fluid  jet  amplifier  is 
shown  with  the  power  jet  attached  to  the  left  boundary  wall. 
In  this  condition,  the  entrainment  flow  through  the  left  con- 
trol port  is  greater  than  that  through  the  right  when  both 
control  ports  are  open  to  the  atmosphere.   If  the  left  con- 
trol port  is  sealed  from  the  atmosphere  while  the  right  port 
is  left  open  and  a  plt)t  of  the  left  control  pressure  and 
flow  is  made,  a  plot  similar  to  Figure  5.2  will  be  obtained. 
This  graph  shows  that  as  the  left  control  pressure  is  in- 
creased above  atmospheric  pressure,  a  point  is  reached  where 
the  jet  is  supplied  all  the  entrainment  flow  it  can  take.   It 
then  becomes  separated  from  the  left  wall  and  switches  to  the 
right.   When  this  occurs,  the  flow  suddenly  decreases  due  to 
a  different  pressure  in  the  interaction  chamber  and  now  pre- 
sented to  the  left  control  port.   If  the  left  control  port 
pressure  is  then  decreased,  a  point  is  reached  where  the  jet 
will  switch  back  to  the  left  side. 

If  the  amplifier  has  perfectly  symetrical  control  ports, 
the  "stream-on-right"  curve  will  represent  the  pressure-f low 
characteristics  at  the  right  control  port  when  the  stream  is 
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attached  to  the  left  wall.  iJ'rom  the  graph  it  is  seen  that  the 
flow  differential  required  for  switching  from  left  to  right  is 
Qa-Qq  and  the  pressure  differential  required  is  T^-j?^ 
Similarly  the  flow  differential  to  switch  from  right  to  left 
is  0,3-0,1,  end  the'  required  pressure  differential  is  ?a-m~^r,# 
The  two  sets  of  flow  and  pressure  differentials  are  not  nec- 
cessarily  equal. 

The  effect  of  loading  on  the  output  legs  of  these  de- 
vices is  generally  to  decrease  the  flow  and  pressures  required 
for  switching. 

5.2  Problem  Statement  and  Objective  of  Experiments 

Because  of  the  nature  of  operation  of  proportional 
amplifiers,  the  useful  signals  from  these  devices  are  elevated 
about  a  bias  pressure  level  greater  than  atmospheric.  •  It  was 
not  clear  what  the  effect  of  this  bias  level  would  be  on  the 
switching  characteristics  of  a  bistable  amplifier  being  driven 
by  a  proportional  amplifier.  This  information  was  necessary 
before  undertaking  the  design  of  the  bistable  reaction  jets 
in  order  that  the  pressure  differential  required  at  the  con- 
trols of  the  proportional  amplifier,  to  effect  switching  of 
the  reaction  jets,  might  be  minimized. 

The  test  procedure  was  to  bias  the  control  ports  of  a 
memory  device*  at  various  pressure  levels  above  atmospheric 
and  then  measure  the  pressures  and  flows  required  to  effect 
switching  at  these  elevated  levels. 


*  A  memory  device  was  selected  for  these  tests  as  this 
is  the  type  of  device  used  for  the  first  stage  of  the  reaction 
jets  discussed  in  the  next  chapter. 
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During  the  tests  the  output  iegat-were  loaded  into  an 
orifice  whose  area  was  approximately  0.75  times  the  receiver 
areas  of  the  device. 


5.3  Test  Results  and  Conclusions 


The  memory  device  selected  for  these  tests  was  a 

I 
Corning  Standard  Size  .jLoad  Insensitive  device*. 

A  graph  of  representative  results  is  shown  in  Figure  5.3, 
It  can  be  seen  that  the  effect  of  the  bias  level  is  to  radi- 
cally increase  the  flow  and  pressure  differential  required  for 
switching.  ' 

.  During  these  tests  it  was  found  that  the  input  charac- 
teristics of  each  control  port  were  quite  different.   There- 
fore, to  more  accurately  show  the  flow  and  pressure  differ- 
entials required  for  switching,  the  pressure-flow  relations 
for  both  legs  are  plotted. 

When  a  bias  level  is  present  at  the  control  ports  of 
a  bistable  amplifier,  the  power  jet  is  being  supplied  more  • 
flow  than  it  can  entrain.  In  this  condition  the  bias  flow 
is  attempting  to  operate  the  device  in  a  proportional  manner. 
However,  due  to  the  close  proximity  of  the  bounding  walls, 
the  jet  will  still  reduce  the  pressure  slightly  on  one  side 
and  attach  itself  to  that  wall.  This  increased  pressure 
differential  required  for  switching  can  not  be  tolerated  in 
this  system  as  it  represents  a  large  hysteresis  in  the 
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switching  of  the  bistable  reaction  jets.  This  problem,  how- 
ever, can  be  eliminated  by  making  the  bistable  amplifier 
large  enough  so  thst  the  entrainment  flow  of  the  power  jet 
is  just  equal  to  the  bias  level  flow. 
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CHAPTER  VI 
POWER  AMPLIFIER  AND  BIST^B^E  REACTION  JETS  ANALYSIS  AND  DESIGN 

6.1  Determination  of  the  Required  Output  Thrust 

So  far  little  mention  has  been  made  of  the  dynamics  of 
the  system  to  be  controlled.  As  mentioned  in  the  introductory 
chapter,  the  system  was  to  approximate  an  ideal  inertia  load. 

The  dynamics  of  such  a  system  are  given  by 

• » 
T    =    ZQ  (6.1) 

where  T  =  the  applied  torque  to  the  system 

J  =  polar  moment  of  inertia 

\j  z  angulsr  displacement 

•  s  differentiation  with  respect  to  time 

If  this  equation  is  integrate^  once  with  respect  to 
time  and  the  initial  angular  rate  is  assumed  to  be  zero,  one 
obtains 

i 

6   =   W   =  ~—  t  (6.2) 

Integrating  equation  (6.2)  with  respect  to  time  and  setting 
the  initial  angular  displacement  to  zero  gives 

0  -    4-  4-        .  (6-3) 


"~3~  sr 
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From  the  design  considerations  discussed  in  the  in- 
troductory  chapter,  it  was  desired  tfcpt  the  system  respond 
fairly  slowly  and  be  able  to  receive  large  angular  displace- 
ments.  On  the  basis  of  these  considerations,  it  was  decided 
that  the  parameter  T/J  should,  be  such  that  if  the  system 
starts  from  rest  under  full  acceleration',  it  must  be  travel- 
ing at  an  angular  velocity  equal  to  1.1  radians  per  second 
at  the  end  of  45  degrees  of  angular  displacement. 

If  equation  (6.3)  is  solved  for  t  and  substituted 
into  equation  16.2)  the  result  upon  rearranging  is 

» 

T    «=    k>2  (6.4) 

Substituting  the  values  of  CO  and Q   into  this  equation  gives 
T/J  =  0.772  sec   .  From  equation  (6.3)  the  time  to  travel 
45°,  with  this  T/J  ratio,  is  1.424  seconds. 

Having  determined  the  T/J  ratio,  an  estimate  of  the 
systems  final  polar  moment  of  inertia  was  made  by  accurately 
determining  the  polar  moment  of  inertia  of  the  rate  sensor  and 
multiplying  the  result  by  a  factor  of  three  to  account  for  the 
remaining  equipment  to  be  added  to  the  system.   The  J  so  cal- 
culated was  found  to  be  7.5  in-ozf-sec  which  results  in  a 
required  torque  of  5.79  in-czf.   If  the  lever  arm  is  assumed 
to  be  5.79  inches,  then  the  thrust  required  is  one  ozf. 

6.2  Introductory;  Theory  of  Staging  Unvented  Bistable 
■amplifiers   ' 

This  system  was  to  develop  the  required  thrust  by 
using  the  change  of  momentum  in  a  reaction  jet.   Since  the 
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power  available  from  the  output  legs  of  the  summer- amplifier 
was  too  small  to  develop  the  required  thrust,  this  signal 
had  to  be  amplified. 

Conpagnuolo^  has  shown  thet  this  power  amplification 
can  be  efficiently  accomplished  by  staging  bistable  amplifiers 
in  ascending  order  of  size  with  each  stage  operating  at  the 
same  supply  pressure.  He  has  pointed  out  that  in  order  to 
make  such  a  cascade,  two  design  principles  must  be  adhered  to: 

1)  For  maximum  pressure  recovery,  the  splitter  dis- 
tances  must  be  kept  small,  preferably  in  the  core  region  or 
the  power  jets  which  in  general  is  approximately  six  nozzle 
widths  down  stream.         : 

2)  In  order  that  the  units  be  capable  of  being  loaded, 
each  stage  must  have  a  high  degree  of  stability  which  implies 
moving  the  splitter  down  stream  to  over  approximately  twelve 
nozzle  widths. 

At  the  present  state  of  the  art,  a  bistable  unit  can 
not  possess  both  the  characteristics  of  high  pressure  re- 
covery (high  efficiency)  and  high  stability.   As  the  splitter 

is  moved  down  stream,  the  pressure  which  can  be  recovered, 

$ 
falls  rapidly.   If  the  splitter  is  moved  up  stream,  then  the 

ability  to  operate  under  load,  (stability*)  falls  off  rapidly. 


*  The  core  region  of  the  power  jet  is  that  region  in 
which  the  velocity  of  the  jet  has  not  been  reduced  from  its 
value  at  the  power  nozzle  due  to  viscous  friction  with  the 
surrounding  fluid. 
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Coapagnuolo  suggests  that  since  the  last  stage  is  the 
most  important  for  maximum  efficiency  that  it  be  mace  with  the 
splitter  located  in  the  core  region  of  the  power-jet.   This 
stage  can  then  be  held  into  the  load  by  making  the  preceeding 
stages  high  stability  devices.   The  result  is  that  the  overall 
cascade  exhibits  high  efficiency  and  stability  under  load. 

In  general  the  instantaneous  flow  gain  of  high  stability 
devices  is  approximately  ten.  This  means  that  it  is  possible 
yto  increase  the  size  of  the  power  nozzle  of  each  stage  by 
roughly  a  factor  of  ten.  It  snould  then  be  possible  to  have 
flow  gains  in  a  three  stage  device  of  at  least  1000  and,  if 
the  efficiency  of  the  cascade  is  good,  power  gains  of  greater 
than  5000. 

6.3  Initial  Experimental  Design 

Figure  6.1  shows  an  output  leg  of  a  bistable  amplifier 
turned  through  a  90  degree  angle.   If  a  control  volume  is 
drawn  as  shown  and  the  momentum  equation  for  steady  state 
operation  is  written  for  the  x-direction,  the  resulting 
thrust  due  to  the  exiting  jet  is  found  to  be: 

T   B   PAeVe2  (6.5) 

Assuming  uniform  velocity  over  the  exit  area,  equation  (6.5) 
can  be  rewritten  as 

T   =  P   ^e  (6.6) 

From  this  last  equation,  knowing  the  thrust  desired  and 
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control  volume. 


I 

Figure  6.1 

Control  Volume  for  Determining  the  Output  Thrust 
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assuming  a  reasonable  value  for  the  exit  area,  one  can  deter- 
mine the  flow  required  to  develop  the  desired  amount  of  thrust. 

If  one  further  knows  the  differential  flow  required 
for  switching  the  first  stage,  then  the  flow  gain  needed  can 
be  readily  determined. 

It  was  shown  in  the  preceeding  chapter  that  it  is  de- 
sirable to  have  the  antra inment  flow  of  the  first  stage 
approximately  equal  to  the  output  bias  flow  of  the  driving 
proportional  amplifier.   This  entra inment  flow  depends  on 

1  2 

the  geometry  of  the  amplifier  and  the  supply  pressure.  It 
can  be  seen  that  if  sufficient  experimental  data  were  avail- 
able that  it  would  be  possible  to 

1)  determine  the  supply  pressure  required  for  a  chosen 
first  stage  amplifier  design  from  the  entra inment  flow  versus 
supply  pressure  characteristics  and 

2)  with  the  supply  pressure  chosen  as  above,  deter- 
mine the  differential  flow  required  for  switching. 

It  would  then  be  a  simple  matter  to  determine  the  flow  gain 
required  and  the  number  of  amplification  stages  needed. 

Lacking  any  such  quantitative  deta,  a  preliminary 
two  stage  bistable  amplifier  was  constructed  in  order  to 

* 

obtain  the  following  information:  (listed  also  are  the  rea- 
sons -this  information  was  desired) 

lj  Entra inment  flows  of  a  first  stage  design  at  vari- 
ous supply  pressures  to  determine  the  supply  pressure  setting 
to  match  the  output  flows  of  the  summer-amplifier. 

2)   Stability  under  lo&d  of  the  unit  to  check  the  set- 
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ting  of  the  splitter  distances,  well  set  back  distences,  etc. 

3)  Output  pressure -flow  relationships  to  determine 
the  optimum  loading  for  maximum  torque  output. 

4)  The  number  of  stages  required  to  develop  the 
desired  torque. 

5)  The  thrust  output  for  various  supply  pressures  and. 

loads. 

The  important  dimensions  of  the  amplifiers  used 

i 
in  this  initial  design  are  as  follows:* 

•  First  ytage 

Power  nozzle  width  *  w  =.  0.020  in. 
Control  nozzle  width  =•  w 
wall  setback    ■»    3w 
Splitter  distance    =.   14w 
Aspect  ratio  »  h/w  =.  1.5w 
Receiver  width  -  3w 

Second  Stage 

.Power  nozzle  width  s  w  »  0.0625  in. 
Control  nozzle  width  *  w 
Wall  setcack   *   1.4w 


*  Many  of  the  dimensions  givyn  here,  auch  as  wall  set' 
back  and  receiver  width,  are  based  on  carefully  considered 
data  given  in  references  12,  13,  and  14.   Time  and  space  lim- 
itations do  not  allow  a  full  discussion  of  the  effects  of 
varying  these  parameters. 
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Splitter  distance     a    6w 

s 

Aspect  ratio  .   -    1.535w  , 

i' 
Receiver  width   s   1.5w 

The  devices  were  n8de  out  of  plexiglass  using  a  pan- 
tograph engraving  w 6 h ine . 

The  results  of  tests  performed  on  this  initial  design 
showed  that  at  a  supply  pressure  of  two  psi  the  entrainment 
flow  in  the  attached  side  was  approximately  0.04  SOFIA.     This 
flow  was  approximately  equal  to  the  no  loaa  output  flow  of 
the  summer-amplifier  when  operating  at  the  supply  pressure 
expected  for  final  design  setting. 

This  device  was  quite  stable  provided  the  output  load 

4  • 

was  relatively  small.   Typical  output  curves  for  the  device 
are  shown  in  Figures  6.2  end  6.3.   In  Figure  6.2  the  device 
is  not  being  driven  by  an  external  signal  and  therefore  can- 
not be  heavily  loaded.   The  overload  area  represents  condi- 
tions under  which  the  second  stage  power  jet  commences  to  oscil- 
late. In  Figure  6.3.  the  device  is  being  driven  by  another 
bistable  amplifier.   It  is  apparent  that  the  device  can  now 
be  more  heavily  loaded  before  the  overload  region  is  reached. 

In  order  to  determine  the  value  of  the  exit  nozzle 
area  for  maximum  thrust,  a  number  of  nozzle  load  curves  were 
drawn  on  a  plot  similar  to  Figure  6.3.   It  was  found  that 
maximum  thrust  output  of  one  nozzle  corresponded  closely  to 
the  point  of  maximum  power  transfer  between  the  amplifier 
cascade  and  the  nozzles.  Since  loading  decreases  the  flow 
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out  of  one  leg,  part  of  the  flow  is  diverted  into  the  other 
leg  where  it  develops  thrust  in  the  opposite  direction.  There- 
fore, considering  the  total  thrust  developed  frcm  the  device, 
it  was  found  that  the  point  of  operation  for  maximum  total 
thrust  ley  slightly  up  the  output  curve  toward  the  region  of 
higher  flows. 

The  results  also  showed  that  in  order  for  this  two 
stage  device  to  develop  an  ounce  of  thrust  that  the  supply 
pressure  would  have  to  be  over  ten  psi.   Since  this  was  in- 
compatible with  the  entrainment  flow  desired,  the  final 
design  would  have  to  have  at  least  three  stages. 

6 .4  Final  Three  Stage  Amplifier  Design 

In  order  to  use  the  -data  gained  from  the  two  stage 
amplifier  to  full  advantage,  it  was  decided  to  use  these 
first  two  stages  in  the  final  design  and  alter  them  only 
slightly  in  order  to  accommodate  the  third  stage.   This 
further  allowed  the  power  nozzle  area  of  the  third  stage  to 
be  calculated  from  the  data  obtained  on  the  two  stage  device. 

The  procedure  used  to  determine  the  dimensions  of  the 
third  stage  are  as  follows:  Assuming  incompressible  flow*, 
the  flow  from  the  third  stage  power  nozzle  is  given  by 


^J3   '=  C<i.  AJ3~\/-§-  Ps       (6-7) 


*  This  assumption  is  reasonably  valid  as  the  final 
device  was  to  operate  at  two  psi. 
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where  0,^3   s     flow  from  third  stage  power  nozzle 
Aj3   =  area  of  third  stage  power  nozzle 
C^    =  discharge  coefficient 
Similarly  the  second  and  first  stage  flows  are  given  by: 


-»« 


4J2   :  °i    A^2-V"P"  **       l6'8) 


Qjx  =     cd  aji-n/"J"  ?s      (6'9)  ' 


The  desired  thrust  is  given  by 


T   =      P°'t  <*jn]  2  (6.10) 


Aj3 

where  it  is  assumed  that  the  load  orifice  areas  are  equal  to 
the  areas  of  the  third  stage  power  jet  and  the  factor  0.8  is  to 
account  for  leakage  flow  out  the  opposite  leg. 

m 

from  the  two  stage  amplifier,  numerical  values  are 
known  for  equations  (6.8)  and  (6.9).   i\irthermore ,  since  Ps 
is  known  and  assuming  a  value  for  c^,  equation  (6.7)  reduces 
to  '• 

<ij3  =   (Constant)  A^3  (6.11) 

If  the  numerical  values  of  equations  (6.8)  and  (6.9) 
are  substituted  into  equation  (6.10)  along  with  equationl6.ll), 
a  quadratic  equation  in  A  13  is  obtained.   Solving  this 
quadratic  equation  gives  the  third  stage  nozzle  area.   Since 
all  dimensions  of  the  third  stage  are  determined  from  the 
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nozzle  width,  assuming  an  aspect  ratio  for  the  third  stage 
will  allow  the  nozzle  width  to  b.e  determined  and  hence  spec- 
ify the  other  dimensions  of* -the  third  stage. 

Following  the  above  procedure  resulted  in  the  following 
dimensions  for  the  final  design: 

First  Stage 

Power  nozzle  width  *  w,  »  0.020  in. 
Control  nozzle  width  -  i.5w^ 
Wall  setback  s   3.0w^ 
Splitter  distance  =   18w^ 
Receiver  width  »   4w^ 

Aspect  ratio  »  1.5W-, 

2 
Nozzle  area  *  0.0006  in. 

Second  Stage 

Power  nozzle  width  3  w2  *  0.0625  in. 

Control  nozzle  width  -  w2 

Wall  setback  •  2.5w2 

Splitter  distance  •  14w2 

Receiver  width  s   3w2 

Aspect  ratio  •   1.535w2 

Nozzle  area   =•  0.00655  in. 

Third  Stage 

Power  nozzle  width  »  w3  »  0.125  in. 
Control  nozzle  width  -  w3 
Wall  setback  »  2wt. 
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Splitter  distance   6w^ 
Receiver  width      1.4w^ 

Aspect  ratio  2.03w3 

p 

Nozzle  area         0,0325  in. 

A  photograph  of  this  device  is  shown  in  figure  D.4. 

15 
The  circular  chambers  are  vortex  impedance  matchers 

modeled  after  those  used  by  Compagnuolo1  . 

In  Figure  6.4  is  shown  the  output  thrust  of  the  ampli- 
fier cascade  versus  supply  pressure  and  input  flov;  rate.   The 
output  flow  of  this  final  design  oscillates  due  to  the  load- 
ing  conditions  unless  the  cascade  is  being  driven  by  an  ex- 
ternal signal. 

Figure  6.5  shows  the  input  impedance  characteristics 
for  the  right  control  port  when  the  first  stage  power  stream 
is  attached  to  the  right  wall.   The  supply  pressure  is  1.96 
psi. 

6 .5  Reaction  Jet  Cascade  Dynamics  and  Transfer  Junction 

The  reaction  jet  cascade  is  characterized  by  two 
phenomena  which  are  important  when  the  device  is  placed  in  the 
control  loop,  of  the  final  system: 

1)  A  hysteresis  loop  is  'present  in  the  switching  of 
the  thrust  from  one  side  to  the  other 

2)  There  is  a  finite  time  delay  between  the  time  a 
switching  signal  is  received  at  the  control  ports  until  the 
thrust  developed  is  switched  from  one  side  to  the  other. 

The  hysteresis  loop  is  due  to  the  characteristics  of 
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bistable  devices  requiring  a  finite  flov;  ana  pressure  dif- 
ferential  to  exist  across  the  control  ports  to  effect  switch- 
ing. The  half  width  of  the  hysteresis  band  was  determined 
from  actual  measurements  with  the  summer-amplifier  set  at  its 
proper  supply  pressure  and  driving  the  reaction  jet  cascade. 
The  half  width  was  found  to  be  0.0685  psi.   The  time  delay 
involves  two  major  phenomena: 

1)  The  time  it  takes  the  fluid  streams  to  traverse 
physical  distances  (transport  time)  end 

2)  The  time  it  takes  to  charge  the  interconnecting 
volumes  to  the  required  switching  levels. 

Actual  measurements  of  this  time  delay  were  not  made. 

From  theoretical  calculations  in  Appendix  B  the  time  delay 

constant  was  found  to  be  0.022   seconds. 

The  transfer  function  of  the  reaction  jet  cascade  is 

v 
thus  given  by 

±  T       Z     e~     3       with  hysteresis  of!  0.0665 
A?  Psi 

(6.12) 
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CHAPTER  VII 
CONTROL  SYSTEM  SYNTHESIS,  PERFORMANCE  AND  ANALYSIS 

7.1  System  Synthesis 

A  major  part  of  the  problem  of  synthesizing  the  system 
components  has  been  accomplished  by  careful  consideration  of 
the  role  each  of  these  components  was  to  play  in  the  final, 
system.  However,  in  order  to  finally  bring  these  components 
together  into  a  working  system  it  is  necessary  to  know  how, 
each  component  will  behave  when  actively  loaded  into  the 
component  it  is  to  drive  in  the  final  design. 

Of  critical  importance  is  the  setting  of  the  power 
supply  pressure  to  the  summer-amplifier.  This  pressure  must 
be  high  enough  to  prevent  the  reaction  jet  cascade  from  oscil- 
lating and,  at  the  same  time,  allow  for  maximum  amplification 
of  the  rate  and  position  control  signals.   It  can  reacily  be 
seen  that  had  the  reaction  jet  cascade  been  designed  with  no 
consideration  for  this  supply  pressure  setting,  that  the  mag- 
nitude of  the  control  signal  required  for  switching  the  jets 
might  have  been  so  large  as  to  preclude  the  operation  of  the 
system. 

The  rate  sensor  had  been  tested  with  its  outputs  loaded 
into  passive  orifices  approximating  the  control  orifices  of  the 
summer-amplifier.  It  was  thus  necessary  to  recalibrate  the 
rate  sensor  with  its  outputs* loaded  into  the  summer-amplifier. 
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For  accurate  results  the  summer-amplifier  must  be  set  at  its 
proper  supply  pressure  and,  ,in  addition,  must  have  the  second 
set  of  controls  activated.   This  is  necessary  as  zhe   input 
impedance  of  the  control  nozzles  is  slightly  affected  by  the 
supply  pressure  and  bias  level  at  the  other  controls.   The 
steady  state  results  of  these  calibration  tests  on  the  rate 
sensor  are  shown  in  Figure  7.1.  For  these  tests,  the  rate 

4 

sensor  output  was  loaded  into  the  summer-amplifier  with  the 
power  supply  pressure  set  at  its  operating  level.   The  second 
set  of  controls  were  not  activated.   The  results  shown  in  this 
figure %were  not  obtained  from  an  accurate  rate  table  and  are 
considered  reasonably  accurate  only  in  the  range  ^0.5  radians 
per  second. 

In  order  to  set  the  maximum  pressure  signal  from  the 
position  sensor*,  it  was  necessary  to  know: 

1)  The  gain  ratio  of  the  two  summer-amplifier  control 

ports  and 

2)  The  final  system  T/J  ratio. 

As  was  pointed  out  in  the  discussion  of  the  signal  comparator, 
the  gain  ratio  depends  on  the  summer-amplifier  power  supply 
pressure  setting  and  control  bias  levels.   Thus,  to  obtain  the 
correct  gain  retio,  a  new  set  of  curves  similar  to  those 
shown  in  Figure  4.1  were  constructed  for  the  proper  supply 
pressure  setting  and  control  bias  levels. 


*  The  maximum  signal  occurs  at  ±90  degrees  of  angular 
rotation. 
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As  shown  in  Appendix  C.l}  the  T/J  ratio  defines  the 

switching  curve  equation.  The  criteria  used  to  set  the  maximum 

pressure  signal  from  the  position  sensor  was  as  follows:   If 

*  » 

the  system  starts  with  an  initial  error  of Qz   90  degrees  and 
Q   *  0  radians  per  second,  the  maximum  rate  signal  at  "Ghe 
point  of  theoretical  optimum  switching  will  occur  at  Q    = 
45  degrees.  To  effect  switching  of  the  jets,  the  position 
signal  should  then  equal  the  rate  signal  at  this  point.  The 

0 

T/J  ratio  defines  the  value  of  this  peak$  •  From  Figure  7.1, 
tha  pressure  signal  corresponding  to  this'value  of  Q    is  deter- 
mined. This  then  is  the  value  that  the  position  signal  must 
have  at  45  degrees.  Knowing  the  shape  of  the  real  switching 
curve  and  the  value  of  the  gain  ratio  of  the  summer-amplifier, 
the  value  of  the  maximum  position  signal  can  readily  be  deter- 
mined.    ' 

In  the  final  system,  a  pair  of  adjustable  orifices 
were  placed  between  the  output  of  the  position  sensor  and 
the  input  to  the  summer- amplifier  to  make  fine  adjustments  of 
the  switching  curve. 

Figure  7.2  shows  a  block  diagram  pf  the  final  control 

system  composed  of  the  various  transfer  functions  from  the 

previous  chapters.  A  summary  of  the  numerical  values  for 

the  parameters  is  given  below: 

"*!  --       1.19   _£si 

rad/sec. 

*KP   s   1.14    P3i 
*  rad 


*These  gains  have  included  the  gains  of  the  summer- 
amplifier.  The  summer  amplifier  multiplied  the  rate  signal 
by  six  and  the  position  signal  by  4.5 
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Figure  7.2 
Block  Diagran  of  Complete  Control  System 
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Figure  7.3  shows  thef  same  system  in  circuit  diagram  form*. 
A  photograph  of  the  complete  control  system  is  shown  in 
Figure  7.4. 

7.2  System  Performance 

To  obtain  performance  data  on  the  system,  a  rotary 
potentiometer  was  connected  to  the  control  system  and  the 
output  signal  displayed  on  an  X-Y  recorder.   The  results  are 
shown  in  Figures  7.5  thru  7.9. 

Figure  7.5  shows  the  response  of  the  system  to  a  num- 
ber of  position  errors.   The  reaction  Jets  switched  in  a  near- 
ly perfect  manner  when  responding  to  counter  clockwise  errors. 
The  response  to  clockwise  errors,  however,  resulted  in  a 
slight  amount  of  "after  end  point  chattering".   The  magni- 
tude of  the  limit  is  seen  to  be  erratic,  averaging  roughly 
8.5  degrees  at  approximately  0.5  cps.   This  erratic  response 
is  due  to  noise  from  the  proportional  amplifiers  in  the 
system.   Considering  the  T/J  ratio  of  this  system,  the 


*  The  symbols  used  in  this  figure  follow  those  recom- 
mended by  3oothe^5. 
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switching 

curve 
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p^sii;ion 
sensor 


Figure  7.3 

System  Circuit  Diagram 
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Figure  7*^- 
Complete  Control  System  and  Demonstration  Modwl 
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respor.se  time  for  the  errors  shown  was  found  to  be  nearly 

0 

optimal. 

Figure  7.6  shows  the  system  operation  with  a  slightly 
lower  summer-amplifier  supply  pressure  setting.   By  lowering 
the  supply  pressure,  the  magnitude  of  the  hysteresis -band  was 
decreased,  resulting  in  the  marked  decrease  in  the  magnitude 
of  the  limit  cycle.   Unfortunately,  the  noise  present  in  the 
system  could  then  switch  the  reaction  jets  more  readily  re- 
sulting in  a  slightly  less  optimal  time  response  to  large 
errors  because  of  random  switching, 

Figure  7.7  shows  the  system  performance  with  the  rate 
sensor  disconnected  from  the  control  loop.  As  pointed  out 
more  clearly  in  appendix  C.2,  the  control  system  should, 
theoretically,  be  completely  unstable.   Because  of  the  finite 
damping  in  the  system,  in  practice  a  new  limit  cycle  was 
reached  at  a  much  higher  amplitide  and  lower  frequency. 

Figures  7.3  and  7.9  show  an  interesting  mode  of  oper- 
ation, possible  with  this  system,   ns  pointed  out  in  the  pre- 
ceding chapter,  the  reaction  jet  cascade  is  unstable  unless 
driven,  due  to  the  heavy  loading  on  th^  output  legs  of  the 
last  stage.   This  back  pressure  is  so  high  that  it  causes  the 
first  and  second  stage  amplifiers  also  to  oscillate.   Intro- 
ducing small  control  flows  into  one  side  of  the  first  stage 
cause  the  pulsing  to  occur  more  often  on  one  side  than  on  the 
other  in  direct  proportion  to  the  magnitude  of  the  input 
signal.   Thus  the  system  effectively  operates  in  a  pulse 


95 


iti;}:::.1:.;, 


r,, 


S     5    }*i 

i  '      i — I — I — H* 


»    —        -        r      .  — ( ► 

?■  ?  3 

»    ,       rf H f 


I  I  _    ■♦.;»_j.j  j     ij 


Pl^urfe  7.7     Syalteo  HteaponBe,  Bfeng-B^ng  Control*  "Wifctkrut  Ratd^eefllnwjk 


IT 


1 


*     £      S      3      R      8 


!    >         ■  t    '    '  |>    '      >i  (  }'     I  I    {     '»     J     I'       'I 
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width  inoculation  mode*. 

In  Figures  7.3  and  7.9  this  pulse  width  modulation  node 
was  obtained  by  lowering  the  summer-amplifier  supply  pressure. 
The  frequency  of  oscillation  is  too  high  (approximately 
100  cps)  for  the  system  to  respond  to  the  individual  torque 
pulses . 

Figure  7.8  shows  the  response  ,t°  various  disturbances 

l 

•.-■■ 

while  operating  in  this  pulsfe  width  modulation  mode.   Of 

particular  note  is  the  disappearance  of  random  fluctuations 

cue  to  noise.  A  certain  amount  of  stiction  is  noted  caused 

by  the  potentiometer  attached  to  the  system. 

Figure  7.9  shows  the  system  response  with  the  rate 

i 
sensor  disconnected  from  the  control  loop,   if  the  load  were 

completely  free  of  camping,  the  response,  of  course,  would 

not  camp  out.   A  compsrison  of  this  figure  with  Figure  7.7 

shows  the  functioning  of  the  vortex  rate  sensor  and  the  part 

it  plays  in  the  stabilization  of  the  system. 

7.3  System  .-.ne lysis 

This   control   system  falls   into   the   category  of   non- 
linear  control  systems   because   of  the   bang-bang  method  of 
controlling  the    load.      The  most    convenient  way  of  analyzing 
this   system,    in  order   to   predict   the  magnitude  and  frequency 
of  the    limit   cycle  and   the    stability  of  the    system,    is   by 
using   the   describing-function  method. 

Gibson  points  out  that  this  approximation  method  is 


Tor  a  more  complete  explanation  of  this  mode  of  oper- 
ation see  reference  17. 
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quite  accurate  in  systems  containing  sufficient  lev/  pass  filter- 
ing.  The  nature  of  an  inertial  load  amply  fulfills  this  re- 
quirement as  it  refuses  to  respond  (filters)  to  high  frequency 
s  i  gna 1 s . 

This  method  essentially  consists  of  assuming  that  yhe 
input  to  the  non-linear  element  is  sinusoidal.   The  resulting 
output  is  then  written  in  the  form  of  ti  Courier  series.   Then 
it  is  assumed  that  only  the  fundamental  component  of  this 
series  is  important  in  relation  to  the  response.   The  coeficient 
of  this  fundamental  frequency,  which  is  dependent  on  the  mag- 
nitude of  input  signal  but  independent  of  the  frequency,  is 
used  as  an  ecuivalent  gain  (K   )  of  the  non-linear  element. 
The  rest  of  the  system  being  linear  is  frequency  dependent 
but  not  dependent  on  the  magnitude  of  the  input  signal.   3y 
plotting  the  equivalent. gain  function  of  the  non-linear 
element  and  the  frequency  response  of  the  linear  elements  on 
a  polar  plot,  a  point  of  intersection  may  be  found  which  re- 
presents a  possible  point  of  operation.   If  the  point  is 
found  to  be  a  stable  point,  then  the  frequency  and  magnitude 
of  oscillation  can  readily  be  determined. 

This  method  is  applied  to  the  analysis  of  this  system 
in  Appendix  (J.  2  The  results  closely  predict  both  the  magni- 
tude and  frequency  of  the  limit  cycle  shown  in  figure  7.5*. 


""The  data  on  which  the  calculations  were  based  was 
obtained  from  the  system  just  after  figure  7.5  was  recorded. 
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APPENDIX  A 


RATE  3ENJCR   CALCULATIONS 


A.l     Rate   Sensor  Design   Calculations 


Given 


A? 


MS 


RQ    =     3  in. 


60 


CO 


max  s 


max 


70°/sgc     i     1.221  rad/sec 
0.1   lbf/in2 


A  $w       ^max    "     3°  =   °-0523  rad« 


2.  34a  10 


-3  lbf-sec2 


ft4 


y 


1.623.;i0    *     ft^/sec 


From  equation   (2.19) 


us^ 


144  A  ? 

*PA# 


-v    /  (144)     (O^lj.     , 
^4(0.00243}  ,(0.0523) 


Us  =   171.2   ft/sec 


jTrom  equation    (2.21) 


*4      " 


irrs     -  4op8  . 


IT  1 0.156 


■h 


62)       -   4(0.0212)     (0.1562), 


J171.2' 


144 


Q,  *   0.151   ft3/sec    9.05  ft3/min. 
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From  equation  (2.22) 


Re  =  -a 


h/2 


KR02 


Re  *   (0-151)  (0.0251)      m  31Q 
(1.623a154)  (0.0625) 

From  reference  3,  page  327,  we  find  that  this  value  of  Re  is 

off  the  scale  but  E  is  leveling  off  to  roughly  0.66.   A 

"factor  of  safety"  was  introduced  here  and  E  was  taken  to 
be  0.5. 

Rearranging  equation  (2.20)  gives 

Ro2  CJ  2 


rpo 

&d  us 

r 

= 

(0.0625)     (1.221)    (0.5) 

PO 

(0.0523)     (171.2) 

r 
po 

s 

4.27  A    10"3ft.           5.12 

-2  . 

i   mc . 

i 


Rearranging  equation  (2.23)  gives 


^\Apo2      "      ^ 


-.      \/(5.12   JL   10-2)2    .iP^0212ll 

2 


a     -      0.0475   in. 
A. 2      Calculation  of   First   otare   oupnly  Pressure 

From  Figure   2.8,    Qc  Qax  z   0.107   in3/sec   and  take   Q,c  ^n*  0. 
Then  V^J^c.  *     §tifl5=T-    -     265  in/sec   , 


lOh 


ifrom  equation    (2.24),        VDs 


P 


wr    t 


»*e 


where              tan/9-     WP^2         -      °-,GG5  ■    =  .05 
u      0.1 


wc     =      0.015  in. 


Wp    s      C.01C   in. 


Therefore,    Yp2        2.43  a  106      in  /sec' 


?rox  equation    (2.23) 


ApVVCcP     AP       V—  PP 


VP2  P/2 


P 


C*    2 


From  Corning  design  data  sheets,  C^p  =x  0.70 


Thus  P  «  0.28  psi  -  7.7  inches  of  water. 
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APPENDIX  B 
CALCULATION  OF   Thi  REACTION  JjET  TIES  LELAY  CONTACT 

The  time  required  for  the  power  jet  to  traverse  the  di-S' 
tance  between  the  nozzle  and  the  receiver  ports  is  given  by 


V 
avg 


where  1  -   splitter  distance 
vavg 


=  the  average  velocity  of  a  particle 
while  traveling  between  the  power 
nozzle  exit  and  the  receiver. 


Reference  19  gives  the  dynamic  pressure  profile  of  a  sub- 
merged  jet  at  various  distances  from  the  power  nozzle.   If 
it  is  assumed  that  the  velocity  decreases  linearly  from  noz- 
zle exit  to  the  receiver,  then  v"av  is  given  by 

vavg        — « £-. 


where     V.        =  ~\/ ~~D~     ^j     =     &°zzle  exit   velocity 

Vs   -   jet  velocity  at  the  entrance  to 
the  splitter 

The  value  of  V3  can  be  determined  from  the  grephs  given  in 
reference  19  if  the  splitter  distance  is  known.   For  the  first 
stage: 

splitter  distance  s  13w  s  0.36  in  *  1 
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V        =  -n    /   ^         Pt       _-v    /laLXUifiJ ^   Z    5   9     '   in3 


P/2    V 


'  in/sec 

18  =      0.45 


P/2   V/ 

V18     c     vs     «      3,96  a  10  3     in/sec 


Vavg        a  V.i  *Vg      =     4.93  a   103     in/s 


ec 


2 


<• 


•       —  1  —         0  .  33 _     „„  .  r\~& 

Similarly  for  the  second  and  third  stages: 

t2  =   197  jC  lb"6  •  sees. 

tr.   -   137  X  10'6  sec3. 

Che  time  required  for  the-  flow  in  the  control  nozzle  of 
an  intermediate  stage  to  build  to  that  flow  level  necessary 
for  switching  can  be  approximated  from  the  time  constant  as- 
sociated  with  the  volume  connecting  the  two  stages  and  the 
resistance  of  the  control  port.   The  time  constant  is  given 
by 

T  =   RC 
where  R  =  impedance  of  the  control  port  lbf -sec/in 

« 

C  =  capacitance  due  to  the  interconnecting 
volumes  in5/l°f 

The  value  of  H  was  determined  by  linearizing  the  input  impedance 
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curves  of  each  nozzle  about  an  assumed  operating  point.   The 
value  of  C  is  given,  vfor  adiabatic  conditions,  by 


C  = 


where  v  =  connecting  volume  (in3) 

p  ■  mass  density  (Ibf-sec^/in**) 

K  =  adiabatic  constant  =  1.4 

R   =  gas  constant  =  2.48  A   105  ( ia2/sec*"0x*) 

T  =  temperature  (°*0 

For  line3  connecting  the  first  and  second  stage: 

R     =   0.237   lbf-sec/in5 

» 

The  value  of  v  was  found  to  equal  0.279  in  •   Therefore 

I 

v      s     (0.279) 

PxrtT        (1.1*9^10?) (1.4)(2.46a10s)(550) 

C   =   0.013   in5/lbf 

Hence 

T  Z        .00372  sees. 

The  relationship  between  the  flow  into  the  connecting  volume 
(assumed  constant)  and  the  flow  through  the  control  nozzle 
(variable)  is  given  by 


^o 


Mn       Ts+1 

The  value  of  %j_n  is  known  from  the  nozzle  geometry  of  the 
first  stage.   ^0  is  zero  before  switching  occurs.  After 
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switching  Q,0  increases  until  it  reaches  the  value  required 
for  switching.  The  value  required  for  switching  was  taken 
to  be  nearly  equal  to  %ia.   Therefore,  the  time  required  for 
the  flow  to  build  to  the  switching  level  is  approximately 
four  time  constants,  or  t-8  0.0140  seconds.  In  a  similar 
manner  the  buildup  time  between  the  second  and  third  stage 
was  found  to  equal  0.006  seconds.   The  final  time  delay  is 
due  to  the  time  required  for  the  output  flow  to  traverse 
the  length  of  the  output  leg.  3y   an  elementary  calculation 


this  was  found  to  equal  636  a  10~°  seconds 


u 


The   total  time   delay   is   equal   to   the   sum  of  all  the 
time   delays   discussed  ebove,   0.022  seconds. 
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APP3NDIX  C 


SYSTSL1  ANALYSIS 


C.l  Deten.iir.ing  the  Equation  of  "ohe  Switching  Curve 

Consider' a  bang-tang  control  system  as  shown  in  figure 
C.l.  The  equation  of  motion  describing  the  operation  of  this 
system  is  given  by 


If 


e 
e 

9* 


Sign$ 


J 


(c.l) 


is   rewritten  as 

dCJ  T     dCO      d(9_  ^     dOJ 
dt        d£      dt  d£) 


and  substituted  for  Q  in  equation  C.l   one  obtains 


J 


CO     d60     _ 

Equation  C.2  C8n  now  be  integrated  resulting  in 


<ign£> 


(C.2) 


6J1 


2 
<x>2 


^L$*Cl  for   (-;  Q  (C.3a) 


rn 


<$>  °2 


for 


C+)  (9 


(,C.3b) 


2         J 
These  curves  represent  two  families  of  parabolas  in  the  phase 

plane  (CJvs  Q  )   where  Q\   and  0%   determine  the  values  of  the 
initial  conditions.  If  these  curves  are  plotted  in  the  phase 
plane,  it  will  be  found  that  the  curves  corresponding  to 
C]_=.  Cg  *  o  pass  through  the  origin.  They  represent  trajectories 
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figure  C.2 
Generalized  Autonomous  Control  System  Containing 
a  Nonlinear  Element 
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such  that,  if  the  system  is  switched  at  the  tine  one  or  the 
other  curve  is  intersected,  a  tine  optinul  reduction  of  errors 
introduced  into  the  system  will  result. 

Thus  the  switching  line  equations  are  gjlven  by 

-  _ *  ■  J     -    seconc  quacrant 

2  switching  curve 

:  +  +  Cj    .  fourth  quadrant 

2'       J  ^  switching  curve 

C.2  Linit  Cvcle  Analysis 

The  block  diagram  shown  in  Figure  7.2  can  be  reduced 
to  the  form  shown  in  Figure  C.2  if  the  following  conditions 
exist: 

1)  The  reference  signal  is  zero 

i 

2)  The  magnitude  of  the  position  signal  stays 
within  the  linear  region  of -the  switching  curve.* 

In  the  figure,  N  represents  the  nonlinear  element  in  the 
system  and  G  is  the  combination  of  ail  the  linear  elements. 

Gibson  shows  that  a  sufficient  condition  fcr  the 
existence  of  sustained  oscillation  is 

•  * 

G(j<Jj  =-_L 


jxeq 


where      Keq  =.  the  equivalent  gain  of  the  non- 
linear element,  which  is  a  func- 
tion of  the  magnitude  of  the  input 
signal  to  the  nonlinear  element. 


This  condition  is  satisfied  as  shown  by  the  results 
in  Figure  7.5. 
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i^eq  for  a  contactor  with  hysteresis  is  given  by  Gibson  as 

Keq  s  g  (K)  +  j  'o(E)  * 


where  E       -    the  magnitude  'of  the   input   to   the 

contactor  represented  by  A^y*-*- A  ?/J 


and 


g(E)  .  JBL    -\/i-  (R_ )2 

b(E)r-4T__       _a_ 
ITS  E 


*, 


where  a      =     the  half  width  of  hysteresis   band 

(psi) 

T    =       the  magnitude  of  the  torque  output 
(in-ozf ) 


For   this   system, 


0.0686  psi 
7.72   in-ozf 


By  assuming  values   for  EsAr^-t-Aiy'    ,    various  values  of  -  l/Keq 
may  be  obtained.      Eor   example 

assume  E  -   3   inches  of  water         Q.108E  psi. 


1 '  (*)  * 

4    (7.72)  -\j\   ^0.0686   \"a 

^    ""      \0.108£/ 


g  = 

Tf(0.1082) 


*  llhese  values  are  normalized  Courier  coeficients. 
The  term  b(E)  is  needed  to  take  care  of  phase  distortion  of 
the  fundamental  harmonic. 
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g   r   70 


'  in-ozf 
psi 


b  -   4Ta  .  m  4(7. 72)(0. 0586 _ 
TTE^      IT  (0.1082)£~ 

Keq  =  89.5  +  -  39.4° 


->-i-  s  0.0112   — E£i 4.  39,4° 

keq  in-ozf 


-57.6 


The  linear  transfer  function  G(s)  is  given  by 


G(S)» 


KjSe"  Tl  V 


"2 


ToS+1 


e 


-T,S 


JS 


2 


in-ozf 
psi 


(C.4) 


Assume  that   the  time  delay  of  the   rate  sensor    (T^)    is 
zero   in  order  to  see  what   effect  rate   feedback  has  on  the 
system.      For  this  case,   the   linear  transfer  function  becomes 


0(3),   *2        (^L32t4l3  4-ljrr3S 


(C.5) 


S2(T2S+1) 


3y  substituting  actual  numerical  values  gives 


0(6)  z   (0.1188) 


r  s2 


J4.89)2    4.89 


2(2.55)  3+1 


e 


-0.0220 


^   (0.043+1) 


(C.6) 
To  show  the  effect  of  rate  feedback  on  the  stability  of  this 
system,  disconnect  the  rate  f-eecback  (K^o)  .  For  this  case 
equation  (0.6)  becomes 

-0.022S 


G(3)  =  0.1188 


S2  (0.046  ■#■!) 


Ilh 


3y  substituting  s*j£j,  the  values  for  the  magnitude  and  phase 

oi'  ^(jCJ;  as  CO   varies  from  zero  to  infinity  ere  obtained.   If 

-his  frequency  response  date  is  plotted  on  a  polar'  plot  along 

with  the  magnitude  and  phase  of  --1 ,  a  plot  similar  to  that 

*-eq 

shown  in  Figure  C.3  will  be  obtained.   It  can  be  seen  that  no 
intersection  of  the  two  plots  occur  r-no  hence  the  system  is 
unstable. 

Substituting  s-jCO   into  equation  (C.oj  and  plotting  this 
function  on  a  polar  plot  (figure  0.4)  shows  the  stabilizing 
effects  of  rate  feedback.   It  can  be  seen  that  a  point  of 
intersection  now  occurs  and  a  stable  limit  cycle  results. 
The  effect  of  rate  feedback  is  to  introduce  phase  lead  into 
the  system  which  "pulls"  the  G(jCl))  curve  away  from  the  -180 
degrees  phase  shift  line  a  sufficient  amount  to  allow  this 
intersection  to  occur.  To  investigate  what  effects  the  time 
delay  in  the  rate  feedback  will  have  on  the  system,  Equation 
(C.4)  may  be  written  in  the  following  form 


G(6)   s  t<L 


o2(Tpo+l) 


If   the   rate   sensor  time   delay  term  is  factored   out   of   the 
quadratic   term  and   combined  with  the   reaction   jet  time   delay, 
the    result   is 

G(o).-    Ko   (  K1T2       o2*-^     o+eri°]        e-^VT3J3 


.2 


(   To   3+lj 

(C.7) 
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Figure  C.3 

Polar  Plot  of  G(jo))  and —  Without  Rate  feedback 

Keq 


Figure  C.*+ 


Polar  Plot  of  a(ico)    and i-   ./ith  Rate  Feedback 

Y 

^eq 
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Comparing  equation   lc.7)   vdth  equation   (C.5J,    the  first  major 

effect  to  be  noted  is  the  increased  transport  tine  delay  given 

by  the  tern  e-^'l* 2^s.   The  other  effect  to  bo  noted  occurs 

in  the  quadratic  term".   If  the  quadratic  terms  given  in  those 

two  equations  are  plotter  over  the  frequency  range  o^GJ^oo, 

it  will  be  found  that  they  differ  in  magnitude'  and  phase  only 

by  a  very  slight  amount  around  the  break  frequency.   Thus,  the 

major  effect  of  the  rate  sensor  time  delay  is  to  increase 

the  transport  time  delay  of  the  over  all  system. 

figure  C. 5  shows  a  Bode  plot  of  equations  (C.5)  and 

(C.7)  where  the  effect  of  the  rate  sensor  time  delay  can  be. 

seen  more  clearly. 

figure  C.6  shows  polar  olots  of  G(jCU)  and . 

r».eq 

The  effect  of  the  rate  sensor  time  delay  is  to  increase  the 
magnitude  of  the  limit  cycle  and  decrease  the  frequency. 
Prom  this  figure,  the  frequency  of  the  limit  cycle  for  the 
system  is  predicted  to  be  9.5  radians  per  second  or  1.51  cps. 
The  predicted  magnitude  can  be  founc  oy   working  back- 
wards through  the  linear  transfer  function  to  get  an  expres- 
sion for  Q   .   This  procedure  gives 


#max  £   r  KT       o 


iZJ.  (j^0)%±i  uon  +  il 


*:2  "2 


where  at  the  point  of  intersection 
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-no 


-  T 


-120 


E  (in.  of  HO) 
CO  (rad/sec) 


ure  C.6 
Polar  Plot  of  G(jco)  and  - 


K 


oq 


119 


CO      =  9.5  radians  per  second. 


Ha king  these  substitutions  and  carrying  out  the  indicated 
calculation  gives 


f)        --   5.68  degrees 


max 
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APPENDIX  D 


PHOTOGRAPHS  OF  COMPONENTS 


Figure  D.1 
Vortex  Hate  Sensor  Showing  Details  of  Signal  Amplification 


^i^-ure  D.2  Summer- Amplifier 
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Figure  D.3 


Angular  Position  Sensor 


(;,.,,,  i 


Figure  D#5f 
Pover  Amplifier  Cascade  and  Reaction  Jets 
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